A111D3  7336b4 


PUBLICATIONS 


NIST  Calibration 
Services  for  Humidity 
IVleasurement 


Contributors: 


A.  Wexler 

R.  W.  Hyland 

S.  Hasegawa 


Coordinator: 


P.  H.  Huang 


U.S.  DEPARTMENT  OF  COMMERCE 
National  Institute  of  Standards 
and  Technology 

Chemical  Science  and  Technology 
Laboratory 

Process  Sensing  Group 
Process  Measurements  Division 
Gaithersburg,  MD  20899 


U.S.  DEPARTMENT  OF  COMMERCE 
Robert  A.  Mosbacher,  Secretary 

NATIONAL  INSTITUTE  OF  STANDARDS 
AND  TECHNOLOGY 
John  W.  Lyons,  Director 


QC- 

100 

U56 


#4677 

1991 


NIST 


NIST  Calibration 
Services  for  Humidity 
Measurement 


Contributors: 

A.  Wexler 

R.  W.  Hyland 

S.  Hasegawa 


Coordinator: 

P.  H.  Huang 

U.S.  DEPARTMENT  OF  COMMERCE 
National  Institute  of  Standards 
and  Technology 

Chemical  Science  and  Technology 
Laboratory 

Process  Sensing  Group 
Process  Measurements  Division 
Gaithersburg,  MD  20899 


October  1991 


U.S.  DEPARTMENT  OF  COMMERCE 
Robert  A.  Mosbacher,  Secretary 

NATIONAL  INSTITUTE  OF  STANDARDS 
AND  TECHNOLOGY 
John  W.  Lyons,  Director 


NIST  HUMIDITY  CALIBRATION  SERVICES 


The  National  Institute  of  Standards  and  Technology  (NIST) , formerly 
National  Bureau  of  Standards  (NBS) , provides  a service  to  Government  agencies 
and  to  the  public  for  the  calibration  of  humidity  measuring  instruments. 
Calibrations  are  performed  by  subjecting  the  instrument  under  test  to 
atmospheres  of  known  moisture  content  produced  by  the  NBS  two -pressure 
humidity  generator.  The  most  accurate  calibrations  are  made  with  the  NBS 
standard  hygrometer,  a device  based  on  the  gravimetric  method.  The  discussion 
that  follows  is  broken  into  two  separate  parts . Part  A describes  the  NBS 
standard  hygrometer.  It  is  based  almost  entirely  on  a 1964  publication  [Ref. 

1] . Pertinent  parts  are  reproduced  here.  Refer  to  [1]  as  indicated  below  for 
more  details,  including  many  references.  Part  B describes  the  NBS  two- 
pressure  humidity  generator. 

Among  the  contributors,  A.  Wexler  and  S.  Hasegawa  are  retired  from  NIST. 
They  have  made  significant  contributions  to  the  original  work  on  the  NBS 
standard  hygrometer  and  the  NBS  two-pressure  humidity  generator,  respectively. 
P.  H.  Huang  has  edited  this  report  by  shortening  the  original  papers  of  A. 
Wexler,  R.  W.  Hyland,  and  S.  Hasegawa  in  order  to  provide  a more  convenient 
reference  for  the  two  original  lengthy  papers. 

A.  NBS  STANDARD  HYGROMETER 

The  NBS  standard  hygrometer  [1]  is  based  on  the  gravimetric  method  of 
moisture  determination.  The  measuring  operation  involves  the  absorption  of 
the  water  vapor  from  a water  vapor- gas  mixture  by  a solid  desiccant  and  the 
determination  of  the  mass  of  this  water  vapor  by  precision  weighing;  it  also 
involves  the  determination  of  the  mass  of  the  associated  gas  of  known  density 
by  counting  the  fillings  of  two  calibrated  stainless  steel  cylinders. 

A sample  of  test  gas  of  constant  moisture  content  from  the  NBS 
two -pressure  humidity  generator  is  allowed  to  flow  through  the  NBS  standard 
hygrometer  and  the  instrument  under  calibration.  At  the  termination  of  a run, 
the  mass  of  water  vapor  removed  by  the  drying  train  is  determined  by  precision 
weighing.  The  drying  train  consists  of  three  interchangeable  glass  absorption 
U-tubes  filled  with  anhydrous  magnesium  perchlorate  and  phosphorous  pentoxide. 
The  total  volume  of  dry  gas  is  measured  by  counting  the  number  of  times  two 
calibrated  cylinders,  each  of  about  .03  m’  capacity,  are  alternately  evacuated 
and  filled  with  the  dry  gas.  By  means  of  a pressure  switch,  a vacuum  pump, 
and  associated  automatic  controls,  each  of  the  two  cylinders  is  evacuated  and 
filled  to  a preset  pressure,  producing  a continuous  flow  of  test  gas.  The 
cylinders  are  immersed  in  a thermostated  oil  bath.  After  each  filling,  the 
gas  collected  in  a cylinder  is  allowed  to  reach  equilibrium  and  then  the  gas 
pressure  and  temperature  are  measured.  The  mass  of  the  dry  gas  is  obtained  by 
multiplying  the  volume  by  the  appropriate  gas  density. 

An  analysis  of  the  random  and  systematic  errors  affecting  the  overall 
accuracy  in  the  determination  of  mixing  ratio  (mass  of  water  vapor/unit  mass 
of  dry  gas)  shows  that  if  0.60  gram(g)  of  water  vapor  is  collected  from  moist 
air,  then  the  estimated  maximum  uncertainty  expected  for  mixing  ratios  between 
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27  mg/g  and  0.19  mg/g  is  12.7  parts  in  10^. 


A-1  COMPUTATIONS.  CORRECTIONS.  SOURCES  OF  ERROR.  AND  ACCURACY/PRECISION 

Mixing  ratio  is  defined  by  the  equation 

r = M/(Vp)^  g (1) 

where 

r = mixing  ratio,  mass  of  water  vapor  per  unit  mass  of  dry  gas,  g/g 

M = mass  of  water  vapor,  g, 

V = volume  of  the  dry  gas  associated  with  the  mass  of  water  vapor,  cm^ 
at  temperature  t and  pressure  B,  and 

p = density  of  the  dry  gas  associated  with  the  mass  of  water  vapor, 
g/cm^  at  the  temperature  t and  pressure  B. 

The  subscripts  t and  B merely  emphasize  the  pressure- temperature 
dependence  of  the  density  and  volume,  and  will  not  appear  in  subsequent 
sections.  The  mixing  ratio  is  a computed  quantity  whose  magnitude  is  the  ratio 
of  two  measurable  quantities  - the  mass  of  water  in  a given  water  vapor- gas 
sample  and  the  associated  mass  of  dry  gas.  The  mass  of  water  is  determined 
directly  by  weighing.  The  mass  of  dry  gas  is  determined  indirectly  through  a 
measurement  of  its  volume  and  a knowledge  of  its  density. 

Each  of  the  measurements  contributing  to  the  value  of  the  mixing  ratio 
as  determined  by  use  of  the  gravimetric  hygrometer  was  considered  in  detail  in 
[1]  and  what  follows  is  a summary  of  that  error  analysis. 

A-1-1  BASIS  AND  NOMENCLATURE  OF  ESTIMATES  OF  ACCURACY 

Although  an  attempt  was  made  to  eliminate  or  reduce  to  a negligible 
value  all  known  sources  of  systematic  error,  there  remain  several  that 
contribute  to  the  uncertainty  in  the  mixing  ratio.  The  nature  and  magnitude 
of  these  will  be  assessed.  All  other  errors  are  assumed  to  be  random  and  are 
treated  as  such.  The  estimate  of  the  accuracy  with  which  the  gravimetric 
hygrometer  can  measure  the  mixing  ratio  of  a gas  sample  is  based  on  both  the 
systematic  and  random  errors.  Standard  deviations  are  computed  and  used  as  a 
measure  of  the  random  errors.  Since  the  mixing  ratio  is  a computed  quantity, 
the  law  of  propagation  of  errors  [2]  is  applied  to  the  random  errors  and  the 
standard  deviation  for  the  mixing  ratio  s(r)  is  given  by  eq.  (2)  (see  [1]  pp. 
397-398).  The  standard  deviation  of  any  parameters,  z,  is  expressed  as  s(z). 

s(r)  - [(|^)2  s2(M)  + 1^)2  s2(V)  + (|^)2  s^Cp)]'/^  . (2) 


Maximum  random  errors  will  be  assigned  as  3 standard  deviations,  and, 
conversely,  if  a maximum  random  error  is  obtained,  its  standard  deviation  will 
be  assigned  as  1/3  of  the  maximum. 

A-1- 2 MEASUREMENT  OF  THE  MASS  OF  WATER  COLLECTED  IN  U- TUBES 
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The  mass  of  the  water  absorbed  by  the  desiccant  in  a U-tube  is  weighed 
against  a tare  tube  of  similar  size,  weight  and  shape  on  an  equal  arm  balance. 
It  depends  not  only  on  the  values  of  the  weights  used  in  the  initial  and  final 
weighings  but  also  on  such  factors  as  the  air  buoyancy  effect  on  the 
U-tube, tare  and  weights;  the  water  vapor  adsorption  on  the  external  surfaces 
of  the  U-tube;  tare  and  weights;  static  charge  on  the  U-tube  and  tare; 
convective  air  currents  within  the  balance  case;  handling  and  treatment  of  the 
U-tube  and  tare;  mass  of  internal  gas  in  the  U-tube;  and  incompleteness  of 
absorption  by  the  desiccant.  The  mass  of  water  vapor,  M,  absorbed  by  the 
desiccant  in  a U-tube  is  given  by 

M = Wf  - W-  + + C3  -h  Cg  + C.  , (3) 


where 

= sum  of  face  values  of  weights  necessary  to  bring  the  balance  to 
equilibrium  after  an  experiment,  with  the  U-tube  on  one  pan  and 
the  tare  on  the  other  pan,  g, 

W-  = sum  of  face  values  of  weights  necessary  to  bring  the  balance  to 

equilibrium  before  an  experiment,  with  the  U-tube  on  one  pan  and 

the  tare  on  the  other  pan,  g, 

Cy  = correction  for  calibration  of  weights,  g, 

C|^  = correction  for  the  effect  of  air  buoyancy  on  U-tube,  tare  and 

weights,  g, 

Cg  = correction  for  water  vapor  absorbed  on  external  surfaces  of 
U-tube,  tare  and  weights,  g, 

Cg  = correction  for  mass  of  gas  sealed  in  the  U-tube,  g,  and 
C-  = correction  for  incompleteness  of  water  vapor  absorption  by  the 
desiccant,  g. 

The  mass  of  water  vapor,  M,  admixed  with  a given  volume  or  mass  of  gas 

is  the  sum  of  the  increases  in  masses  of  the  first  two  tubes  of  the  three -tube 

main  absorption  train,  that  is 


M = m.|  + m2 


(4) 


where 

m.j  = increase  in  mass  of  the  first  tube,  g,  and 

m2  = increase  in  mass  of  the  second  tube,  g. 

The  third  tube  is  a guard  against  moisture  entry  from  back  streaming 

sources . 

About  0.6  gram  of  water  is  collected  in  the  first  tube  and  0.2  mg  is 
collected  in  the  second  tube  in  a typical  mixing  ratio  determination. 

Each  of  the  correction  terms  C^,  Cj^,  and  Cg  given  in  eq.  (3)  is  the 
difference  between  the  corresponding  correction  for  the  initial  and  final 
weighings;  for  example,  where  is  a buoyancy  correction 

applied  to  the  initial  weighing  and  is  a buoyancy  correction  applied  to 
the  final  weighing. 
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The  sum  of  the  corrections  applied  to  the  difference  - W-  for  any  tube 
is  of  the  order  of  2 mg.  The  corrections  constitute  nominally  one -third  of 
one  percent  of  the  total  weight. 

The  standard  deviation  of  M,  s(M),  is 

s(M)  = [s2(m^)  + s2(m2)]^/^  . (5) 


The  above  corrections,  the  errors  involved  in  applying  these  corrections , and 
uncertainties  due  to  other  factors  are  now  discussed  in  detail  (see  [1]  pp. 
399-413) . 

A-1-2-1  Balance  and  Weights  (see  FH  p.399).  The  100 -gram  capacity  equal -arm 
semi-micro  balance  used  to  weigh  the  U-tubes  has  a reciprocal  sensitivity  of 
approximately  0.02  mg/division.  It  is  possible  to  estimate  to  0.1  division, 
that  is,  to  the  nearest  0.002  mg.  Pointer  deflections  are  observed  through  a 
telescope.  The  length  ratio  of  the  right  to  left  balance  arm  is  1.000003. 
Single  transposition  weighings  are  used,  which  compensate  for  the  inequality 
in  armlength.  Based  on  repeated  weighing  experiments  with  stainless  steel 
weights  of  50  and  100  grams,  it  is  estimated  that  the  variability  of  the 
balance,  that  is,  the  standard  deviation  of  a single  transposition  weighing, 
is  s(W^)  = s(W-)  = 0.011  mg. 

A set  of  class  M weights  [3]  is  used  with  the  balance.  The  weights  of 
face  values  1 to  100  grams  were  fabricated  from  Brunton  metal  and  have  a 
stated  density  at  20®C  of  7.89  g/cm^  . The  weights  of  face  values  1 mg  to  500 
mg  were  fabricated  from  an  alloy  of  80  percent  Ni  and  20  percent  Cr 
composition  and  have  stated  density  at  20° C of  8.39  g/cm^.  The  maximum 
uncertainty  in  density  is  estimated  to  be  0.03  g/cm^.  This  contributes  a 
systematic  but  negligible  error. 

The  standard  deviation  of  the  correction  for  the  calibration  of  the 
weights  for  the  single  weighing  is  nominally  5.3  x 10"^  g,  and  for  a 
differential  weighing  s(Cy)  is  ^ 2 (5.3  x 10'^)  g or  7.5  x 10'^  g. 

A- 1-2 -2  Buoyancy  Correction  (see  [11  pp . 399-406).  With  an  equal -arm  balance, 
the  use  of  a tare  is  normally  the  preferred  procedure  for  precision 
differential  weighing  since  it  tends  to  compensate  not  only  for  the  effect  of 
air  buoyancy  on  the  U-tube  and  tare,  but  also  for  other  possible  changes  in 
tube  weight  due  to  handling,  dust  accumulation,  and  moisture  adsorption  on  the 
external  surfaces.  However,  since  the  external  volumes  of  the  tare  and  U-tube 
are  not  precisely  equal,  and  since  barometric  pressure  changes  may  on  occasion 
introduce  significant  density  changes,  a buoyancy  correction  is  applied  to  the 
differential  weighing.  This  correction  is  given  by 
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where 


p.  = density  of  the  ambient  air  at  the  time  of  the  initial  weighing, 
g/cm^ , 

= density  of  the  ambient  air  at  the  time  of  the  final  weighing, 
g/cm^ , 

Vj  = external  volume  of  the  tare,  cm^, 

= external  volume  of  the  unknown  U-tube,  cm^, 

W|^.  = calibrated  value  of  weights  1 g and  above  at  the  time  of  the 
initial  weighing,  g, 

= calibrated  values  of  weights  1 g and  above  at  the  time  of  the 
final  weighing,  g, 

Wg.  = calibrated  value  of  weights  less  than  1 g at  the  time  at  the  time 
of  the  initial  weighing,  g, 

Wg^  = calibrated  values  of  weights  less  than  1 g at  the  time  of  the 
final  weighing,  g, 

= density  of  weights  1 g and  above,  g/cm^,  and 
Dg  = density  of  weights  less  than  1 g,  g/cm^. 


Although  C|^  can  be  as  large  as  0.5  mg,  it  normally  is  of  the  order  of  0.15 
mg.  Each  of  the  terms  in  eq.  (6)  is  summarized  in  the  following  subsections. 


A-1-2-2-1  Air  Density  in  the  Balance  Case  (see  [11  no.  400-402).  The  density 
of  the  air  in  g/cm^  is  computed  from  the  relation 


273.16  (B-0. 003780  e,  RH) 

P = Po  X X J ^ 

T 760 


(7) 


where  T is  the  absolute  temperature  in  K,  B is  the  barometric  pressure  in 
Torr,  e^  is  the  saturation  vapor  pressure  of  water  at  T in  mm  Hg  at  standard 
gravity  and  0°C,  RH  is  the  relative  humidity  in  percent  at  T,  and  p^  is  the 
density  of  air  at  a pressure  of  760  mm  Hg  and  a temperature  of  273. 16K,  and 
has  a numerical  value  of  1.29304  x 10’^  g/cm^  (A-1-4-1) . 

The  standard  deviation  of  the  air  density  is  0.32  x 10'^  g/cm^  and  was 
computed  at  nominal  values  T = 298K,  B = 750  mm  Hg,  RH  = 40  percent,  e^  = 23.8 
Torr,  and  p = 1.167  x 10'^  g/cm^.  The  steps  are  listed  in  table  1 together 
with  the  result. 

A-1-2-2-2  External  Volume  of  U-Tube  (see  [11  pp . 402-405).  The  external 
volume  of  each  U-tube  was  determined  by  weighing  in  air  and  in  distilled 
water.  An  equal -arm  balance  was  mounted  on  a table,  over  an  open  cylinder  of 
distilled  water.  A stainless  steel  wire  was  attached  to  the  right  pan  and  a 
29.4-gram  weight  was  suspended  from  the  end  of  the  wire  and  immersed  in  the 
water.  A U-tube  was  then  placed  on  the  right  pan  and  counter -balanced  with 
weights  on  the  left  pan.  This  provided  the  U-tube  weight  in  air.  The  U-tube 
was  then  attached  to  the  wire,  immersed  in  the  cylinder  of  water,  and  again 
weighed.  During  these  operations,  the  glass  stopcocks  were  kept  in  the  closed 
position.  Precautions  were  taken  to  insure  that  the  U-tube  side  arms  were 
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filled  with  water  and  did  not  contain  trapped  air. 

Consider  the  case  where  the  U-tube  is  weighed  in  air.  Let  be  the 
weight  suspended  from  the  left  pan  and  the  weight  suspended  from  the  right 
pan  via  the  wire  and  immersed  in  the  water.  The  weight  and  the  U-tube  are 
buoyed  up  by  displaced  air  of  density  while  the  weight  is  buoyed  up  by 
displaced  water  of  density  d^ . If  the  density  of  is  given  by  D and  of 
by  Dq,  and  if  the  external  U-tube  volume  at  the  ambient  temperature  is  given 
by  V^,  then  it  can  be  shown  that 


V 


o 


- W2  W2  (P2/D)  - (p^/D) 

+ 

(d2'Pl)  ^^2  ■ Pi) 


(yPp)  (^2  - d^) 
(^2  ■ Pi) 


(8) 


where  subscript  2 refers  to  the  final  conditions,  i.e.,  the  weighing  operation 
with  the  U-tube  suspended  in  water. 

Nominal  values  W.j  = 98  grams,  W2  = 33  grams,  = 29.4  grams,  d^  = d2  = 
0.997  g/cm^,  D = 8.4  g/cm^,  = 11.3  g/cm^.  Pi  = P2  = 0.00116  g/cm^,  d2  - d^  = 
0.00025  g/cm^,  and  the  standard  deviations  for  each  independent  variable  yield 
0.017  cm^  as  the  computed  standard  deviation  of  the  external  volume  s(Vq). 

The  computations  and  result  are  shown  in  table  2. 

A-1-2-2-3  Error  in  the  Buoyancy  Correction  (see  [11  pp.  405-406).  An 
estimate  may  now  be  made  of  the  standard  deviation  of  the  buoyancy  correction 
s(C|^).  For  this  computation,  the  following  nominal  values  of  the  independent 
variables  will  be  used:  ^ Pf  Pi  0.00116  g/cm^,  = W|_.= 

10  grams,  = 7.9  g/cm^,  Dg  = 8.4  g/cm^,  p^  - p-  = 10"^  g/cm^,  and  Wg^  = Wg.= 

0 . 5 grams . 

The  calculations  and  the  result  are  shown  in  table  3 . The  experimental 
standard  deviation  of  the  U-tube  external  volume,  0.04  cm^,  is  used  in  the 
calculations,  since  it  is  the  larger  of  the  calculated  and  experimental 
values.  The  standard  deviation  of  the  buoyancy  correction  is  7.6  x 10’^  gram. 
A-1-2-3  Water  Vapor  Adsorption  on  External  U-tube  Surface  (see  fll  p.  406). 
Water  vapor  is  adsorbed  by  metal  and  glass  surfaces  in  proportion  to  the 
ambient  relative  humidity.  Since  the  balance  is  symmetrical  in  its 
construction,  the  adsorption  on  the  arms,  pans  and  other  components  will  tend 
to  compensate  and  may  therefore  be  neglected.  The  adsorption  on  the  metal 
weights  is  also  negligible  in  magnitude.  At  50  percent  RH  and  nominally  room 
temperature,  "Pyrex"  is  assumed  to  absorb  0.11  pg/cm^,  with  a possible 
systematic  error  of  ±0.11  pg/cm^.  If  the  correction  were  ignored  entirely,  a 
maximum  possible  systematic  error  of  ±0.2  pg/cm^  would  be  introduced.  Since 
the  estimated  maximum  difference  between  external  U-tube  area  and  tare  tube 
area  is  20  cm^,  and  since  the  laboratory  relative  humidity  is  kept  below  50 
percent,  it  is  seen  that  the  correction  may  be  ignored,  leading  to  a maximum 
systematic  error  in  U-tube  weight  of  ±4  pg  due  to  water  vapor  adsorption  on 
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the  external  surfaces . 


A-1-2-4  Static  Charge  on  U-Tube . Static  charge  can  accumulate  on  the  pans  of 
the  balance  and  on  the  U-tube.  The  resultant  force  that  is  exerted  on  the  arm 
of  the  balance  will  introduce  an  error.  To  avoid  or  reduce  any  accumulation 
of  charge,  several  precautionary  measures  are  observed.  The  metal  balance 
case  is  grounded  so  that  charge  on  the  pan,  or  other  metal  parts,  will  tend  to 
leak  off.  The  glass  panels  and  doors  of  the  case  are  metalized.  Strontium  90 
is  kept  within  the  balance  case,  emitting  beta  rays  which  ionize  the  air  in 
the  immediate  vicinity  of  the  U-tube.  In  the  presence  of  ionized  air  there 
is  a greater  tendency  for  any  charge  on  the  U-tube  to  leak  off.  Since  static 
charge  is  evidenced  by  irregularities  in  the  pointer  swing,  no  reading  is  made 
until  the  swing  is  regular.  There  may  still  be  residual  static  charge;  the 
error  this  contributes,  if  any,  cannot  be  readily  estimated. 

A-1-2-5  Handling . Whenever  a U-tube  is  manipulated,  as  in  opening  or  closing 
its  stopcocks,  or  in  inserting  it  in  the  drying  train,  a clean  chamois  skin  or 
lens  tissue  is  used  to  avoid  bringing  the  fingers  and  hands  of  the  operator 
into  direct  contact  with  the  surfaces  of  the  tube.  After  removal  from  the 
drying  train, the  side  arms  are  cleaned  with  chamois  skin  or  tissue  to  remove 
any  residue  from  the  "0”  rings.  The  U-tube,  when  not  in  the  drying  train,  or 
in  the  balance  case,  is  kept  in  a covered  plastic  box  to  reduce  any  tendency 
for  dust  or  dirt  to  accumulate  on  the  external  surfaces.  It  is  dusted  with  a 
camel's  hair  brush  prior  to  suspension  from  one  arm  of  the  balance.  The  tare 
is  treated  similarly.  The  magnitude  of  the  error  arising  from  handling  is 
difficult  to  estimate  and  no  attempt  will  be  made  to  do  so  here,  although  it 
may  be  inferred  from  the  experimental  data  subsequently  presented  (section  A- 
1-2-10)  that  handling  is  one  the  largest  factors  contributing  to  the  weighing 
error . 

A- 1-2- 6 Mass  of  Internal  Gas  in  U-Tube  (see  [11  p.  407).  The  internal  volume 
of  an  empty  U-tube  is  about  50  cm^.  When  filled  with  desiccant,  the 
unoccupied  volume  is  about  30  cm^.  If  air  occupies  this  volume, then  the 
weight  of  this  air  is  part  of  the  weight  of  the  sealed  U-tube.  Not  only  can 
the  temperature  and  pressure  of  the  air  at  the  time  the  U-tube  is  sealed 
differ  for  the  initial  and  final  weighings,  but  the  internal  unoccupied  volume 
will  decrease  due  to  the  absorption  of  water  vapor  by  the  desiccant.  Unless  a 
suitable  correction  is  applied,  an  uncertainty  in  the  difference  between 
initial  and  final  U-tube  weights  is  thereby  introduced.  This  correction  may 
be  as  large  as  1.6  mg. 

The  correction  that  must  be  applied  to  a differential  U-tube  weighing 
due  to  a changes  in  internal  gas  density  and  desiccant  is 

Cj  = V,p,  - V,Pi  - Vj(p,-p,)  - vpf  , (9) 


where 

^ net  internal  U-tube  volume  at  the  time  of  the  final  weighing,  cm^, 

V-  = net  internal  U-tube  volume  at  the  time  of  the  initial  weighing,  cm^, 

= density  of  the  gas  (usually  air)  in  the  U-tube  at  the  time  the  U-tube  is 
sealed  in  preparation  for  the  final  weighing,  g/cm^. 
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p.  ^ density  of  the  gas  (usually  air)  in  the  U-tube  at  the  time  the  U-tube, 
is  sealed  in  preparation  for  the  initial  weighing,  g/cm^,  and 
V decrease  in  internal  volume  due  to  absorbed  water,  cm^. 

A-1-2-6-1  Density  of  Gas  in  U-Tube . The  density  of  the  gas  (taken  here  as 
dry  air)  sealed  in  the  U-tube  is  computed  from  eq.  (7)  in  which  the  relative 
humidity  is  zero.  The  standard  deviation  of  the  air  density,  estimated  in  a 
similar  fashion  to  that  described  in  A-l-2-2~l  is  0.29  x 10"^  g/cm^. 

A-1-2-6-2  Initial  Internal  U-Tube  Volume  (see  [11  dp.  407-408).  To  correct 
for  the  weight  of  air  sealed  in  a U-tube,  it  is  necessary  to  know  the  volume 
occupied  by  the  air.  If  the  U-tube  is  weighed  separately  with  air  and  then 
with  hydrogen,  the  initial  internal  volume  can  be  computed  using  the  equation 


Pa  - Ph 


(10) 


where 

Wg  = mass  of  U-tube  filled  with  air,  g, 

= mass  of  U-tube  filled  with  hydrogen,  g, 

Pg  = density  of  air  in  U-tube,  g/cm^,  and 
Pj^  ^ density  of  hydrogen  in  U-tube,  g/cm^. 

The  standard  deviation  of  the  internal  air  density  is  0.29  x 10’^  g/cm^. 

An  analysis  for  the  internal  hydrogen  density,  the  steps  of  which  are  shown  in 
table  4,  yields  a standard  deviation  of  2 x 10’®  g/cm^. 

For  the  purpose  of  this  computation,  the  experimentally- determined 
standard  deviation  of  the  mass  change  of  a single  U-tube  based  on  flushing 
with  dry  gas,  given  in  section  A- 1-2 -10  will  be  used.  Thus  s(Wg)  =-  "" 

0.08  mg. 

The  standard  deviation  of  the  initial  internal  U-tube  volume  is 
s(V.)  = 0.11  cm^,  as  shown  in  table  5. 

A-1-2-6-3  Change  in  Internal  U-Tube  Volume.  The  internal  volume  of  a U-tube 
is  decreased  between  initial  and  final  weighings  because  of  the  adsorbed 
water . 


The  first  U-tube  of  the  drying  train  absorbs  on  the  order  of  0.6  grams 
of  water  vapor.  This  tube  contains  magnesium  perchlorate.  When  a gram-mole 
of  Mg(C{0^)2  absorbs  water,  the  Mg(CC0^)2  may  combine  with  2,  4 or  6 gram-moles 
of  water  [17]  depending  on  the  water  available. 

The  reaction  is 

Mg(C(!0^)2  + nH20  > Mg(C^0^)2  • nH20  (11) 


8 


where  n = 2,  4,  or  6 . 

Under  the  dynamic  conditions  existing  in  the  case  of  a flow  of  moist  gas 
through  the  desiccant,  it  is  possible  for  the  reaction  to  yield  all  three 
hydrate  forms . 

The  internal  volume  decrease  is  given  by 
V = x(W^  - W.)  , (12) 

where  x is  0.819  cm^/g,  is  the  final  U-tube  weight  in  grams  and  W-  is  the 
initial  U-tube  weight  in  grams.  The  magnitude  of  v is  0.5  cm^.  The  value  of 
0.819  cm^/g  is  obtained  by  assuming  the  true  volume  change  per  gram  of  water 
absorbed  is  midway  between  the  value  for  the  hexadydrate  and  that  for  the 
dihydrate . 

The  amount  of  water  absorbed  by  the  second  and  third  tubes,  which 
contain  P205>  is  less  than  0.2  mg.  It  can  be  shown  that  the  volume  change 
from  the  small  amount  of  water  absorbed  is  negligible  [1]. 

A-1-2-6-4  Error  in  Internal  Gas  Mass.  The  random  uncertainty  in  Cg  may  now 
be  estimated.  The  details  of  the  computations  are  given  in  table  6.  Nominal 
values  of  30  cm^  for  and  V.  and  of  1.17  x 10'^  g/cm^  for  and  p-  were  used. 
The  standard  deviation  s(Cg)  so  obtained  is  7.9  x 10’^  gram. 

A-1-2-7  Incompleteness  of  Water  Vapor  Adsorption.  Bower  [18]  has  shown  that 
the  residual  water  vapor  in  a gas  stream  emerging  for  a U-tube  filled  with 
anhydrous  Mg(C^0^)2  is  0.2  x 10’^  at  a flow  rate  of  1 to  5 liter  per  hour, 

while  Morley  [19]  has  similarly  demonstrated  that  an  absorption  tube  packed 
with  anhydrous  P2O5  will  remove  all  but  0.25  x 10*^  mg  of  water  vapor  per  liter 
of  air  from  a gas  stream  at  a flow  rate  of  three  liters  per  hour.  The  flow 
rate  used  in  the  gravimetric  hygrometer  may  be  as  large  as  2 liter  per  minute. 
Because  of  the  large  difference  in  flow  between  Motley's  work  and  the  present 
work,  no  correction  will  be  made,  and  the  systematic  error  in  the  weight  of 
water  arising  from  the  incompleteness  of  absorption  will  be  conservatively 
assigned  a maximum  value  of  0.5  x 10'^  nig/C  of  air  collected, 

A-1-2-8  Random  Errors  in  the  Measurement  of  Mass  of  Water  Vapor.  A numerical 
estimate  will  now  be  made  for  the  random  error  involved  in  the  measurement  of 
the  mass  of  water  collected  in  a U-tube  based  on  the  uncertainties  of  the 
various  independent  parameters  discussed  above.  The  standard  deviation  s(m) 
is  given  by  eq.  (5).  The  individual  sources  of  random  error,  the  magnitudes 
of  their  contributions,  and  the  computed  standard  deviation  s(m)  are 
summarized  in  table  7, 

The  standard  deviation  of  the  mass  change  of  the  first  tube  is  81  x 10’^ 
gram,  and  of  the  second  and  third  tubes  is  25  x 10'^  gram.  The  figures 
represent  the  error  in  mass  gain  in  a tube  where  both  the  initial  and  final 
masses  are  each  based  on  a single  weighing. 

A-1-2-9  Systemic  Errors  in  the  Mass  of  Water  Vapor.  There  are  two  non- 
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negligible  systematic  errors  which  must  be  considered.  The  first,  arising 
from  the  differences  in  surface  area  of  the  U-tube  and  tare  and  the  water 
vapor  absorbed  on  these  surfaces,  was  discussed  in  A-1-2-2-4  and  A-1-2-3.  The 
maximum  error  expected  from  this  effect  was  ±4  pg  per  tube,  and  since  two 
tubes  are  used  for  the  water  vapor  mass  determination,  the  maximum  error  for 
the  water  vapor  mass  determination  is  ±8  pg. 

The  other  systematic  error  was  discussed  in  A-1-2-6-3.  The  maximum 
uncertainty  in  the  constant  x of  eq.  (12),  important  only  for  the  first  tube 
in  the  drying  train,  was  ±0.057  cm^/g.  Assuming  that  0.6  g of  water  is 
absorbed  in  this  first  tube,  it  is  seen  that  a maximum  uncertainty  of  ±0.034 
cm^  arises  in  the  internal  volume.  Since  the  nominal  air  density  in  the  tube 
is  1.17  mg/cm^,  the  systematic  error  in  the  mass  of  internal  air  and  hence  in 
the  mass  of  water  vapor  collected  may  be  as  much  as  ±0.04  mg. 

The  total  possible  systematic  error  in  the  mass  of  water  vapor  is 
therefore  0.04  + 0.008  = ±0.048  mg. 

A- 1-2 -10  Blank  Runs . Two  types  of  experiments  were  performed  to  provide  a 
check  on  the  computed  estimate  of  the  error  in  the  determination  of  the  mass 
of  water  vapor  absorbed  by  the  main  drying  train.  These  experiments  traced 
the  degradation  in  accuracy  of  the  weighing  process  with  increasing  complexity 
of  the  handling  and  manipulating  operations. 

The  first  experiment  consisted  of  blank  runs  in  which  the  actual 
conditions  of  a mixing  ratio  determination  were  simulated,  except  that  air 
which  had  been  successively  predried  by  Mg(CJ0^)2  and  P2O5  passed  through  the 
U-tubes.  Each  desiccant- filled  U-tube  in  this  experiment  was  weighed,  removed 
from  the  balance,  stored,  and  subsequently  inserted  into  the  main  drying  train 
along  with  one  or  two  other  filled  U-tubes.  Its  stopcocks  were  opened,  it  was 
flushed  with  the  predried  air  at  a flow  rate  of  about  1.7  Jpm,  its  stopcocks 
were  closed,  and  it  was  reweighed. 

The  tube  was  weighed  at  least  twice  (generally  three  times  or  more) 
prior  to  a run,  from  which  an  average  initial  weight  was  computed.  An 
identical  procedure  yielded  an  average  final  weight.  The  average  final  weight 
of  one  run  was  used  as  the  initial  average  weight  of  the  following  run.  Since 
the  change  in  mass  should  be  zero,  the  difference  between  the  initial  and 
final  weight  is  a measure  of  the  error.  Three  different  U-tubes  were  used  to 
form  a total  of  fifteen  runs.  The  results  are  shown  in  table  8. 

The  losses  and  gains  in  weight  (differences)  are  almost  equal  in  number 
and  the  sum  of  these  is  essentially  zero.  The  blank  runs  did  not  disclose  any 
systematic  error.  After  applying  a method  for  pooling  data  described  by 
Youden  [12] , and  assuming  that  these  results  may  be  extrapolated  to  actual 
mixing  ratio  determinations,  it  is  estimated  that  the  gain  in  mass  in  a single 
U-tube  can  be  measured  with  a standard  deviation  of  0.080  mg  and  in  two  tubes, 
0.113  mg. 

A second  experiment  was  performed  in  which  a given  mass  of  water  was 
added  to  a stream  of  pre- dried  air  and  was  subsequently  removed  by  the  drying 
train.  The  water  source  was  a U-tube  containing  about  1 gram  of  distilled 
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water  through  which  the  air  passes  at  a rate  of  flow  of  approximately  1.5  Cpm. 
The  mass  of  water  lost  by  the  source  was  compared  to  the  gain  in  mass  by  the 
first  two  tubes  in  the  drying  train.  The  third  tube,  in  accordance  with  the 
practice  adopted  in  this  work,  was  used  only  as  a guard  tube. 

Each  of  the  U- tubes,  the  one  containing  water,  and  the  three  in  the 
drying  train  containing  desiccant  were  weighed  from  two  to  five  times  before 
and  after  a run  to  provide  average  values  for  the  initial  and  final  weights. 

Five  runs  were  made.  The  results  are  given  in  table  9.  The  differences 
between  the  mass  increase  in  the  drying  train  and  the  mass  decrease  in  the 
water  source  show  roughly  equal  numbers  of  plus  and  minus  signs,  indicating  no 
pronounced  tendency  that  could  be  ascribed  to  systematic  effects.  As  in  the 
previous  experiment,  the  differences  will  be  assumed  due  to  random  errors. 

The  standard  deviation  of  a single  determination  of  the  difference  between  the 
water  available  and  that  gained  is  0.16  mg.  This  is  consistent  with  the 
results  of  the  dry  gas  experiment  which  yielded  0.080  mg  as  the  standard 
deviation  for  the  mass  determination  of  a single  U-tube.  Since  three  tubes 
were  involved  in  each  of  the  present  runs  (one  tube  with  water  and  two  tubes 
with  desiccant),  the  predicted  standard  deviation  of  a run,  based  on  0.080  mg 
is  Vs  X 0.080  mg  or  0.14  mg.  Inversely,  it  could  be  stated  that  the  final 
experiment  yields  a standard  deviation  of  0.092  mg  for  the  mass  determination 
by  a single  U-tube  in  the  main  drying  train  so  that  the  standard  deviation  of 
the  sum  of  the  water  vapor  mass  determinations  by  the  first  and  second  tubes 
is  V’2 (0.092)  =0.13  mg.  The  standard  deviation  of  the  water  vapor  mass 
determination  for  any  run  therefore  will  be  assumed  to  be  0.13  mg. 

A-1-3  MEASUREMENT  OF  THE  GAS  VOLUME  (see  [1]  pp. 413-418) 

The  internal  volumes  of  the  two  cylinders  were  measured  by  weighing  each 
cylinder  empty  and  then  filled  with  distilled  water.  Each  cylinder  volume  was 
computed  using  the  equation 


where 

M = mass  of  water  required  to  fill  the  cylinder  at  temperature  t^, 
grams , 

d = density  of  water  at  temperature  t^  and  pressure  B,  g/cm^, 

= volume  at  temperature  t^  and  pressure  B,  cm^,  and 
= correction  for  water  trapped  in  cylinder  valves  during  filling 
operation,  cm^. 

Five  independent  determinations  were  made  of  each  cylinder  volume. 

Since  the  ambient  conditions  were  not  controlled,  each  determination  involved 
a different  temperature  and  barometric  pressure.  At  atmospheric  pressure,  the 
change  in  cylinder  volume  due  to  barometric  fluctuations  is  negligible;  the 
change  because  of  temperature  fluctuations,  on  the  other  hand,  is  significant. 
Each  experimental  value  of  cylinder  volume  was  therefore  adjusted  to  the 
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volume  it  would  have  at  25° C,  and  these  adjusted  volumes  were  averaged  to  give 
a mean  value . 

Thus 


= V + 

OCX 


where 


(14) 


= volume  at  25 °C  and  atmospheric  pressure,  cm^,  and 
= correction  for  adjusting  each  experimentally  determined  volume  to 
that  of  25°C. 


The  mean  cylinder  volume  is 


V. 


n 


(15) 


where 


EVq  = sum  of  the  experimentally  determined  and  adjusted  volumes,  cm^. 


and 

n = number  of  experimental  determinations . 

The  calibration  results  are  given  in  table  10.  Each  cylinder  has  a 
volume  capacity  of  about  29.7  liters  at  25°C. 

In  an  actual  run  in  which  a moisture  determination  is  being  made,  the 
cylinders  may  not  be  at  the  standard  temperature  of  25° C.  It  will  be 
necessary  to  adjust  the  mean  cylinder  volume  to  that  of  the  test  temperature 
by  applying  a correction  whose  magnitude  will  depend  on  the  difference  between 
these  two  temperatures,  thus 

V = \ + (16) 

where 

V = volume  of  the  cylinder  test  temperature,  cm^,  and 

= correction  for  adjusting  the  mean  cylinder  volume  to  that  at  the 
test  temperature,  cm^. 

_ An  analysis  will  now  be  made  to  obtain  an  estimate  of  the  uncertainty  in 

Vq  based  on  all  known  sources  of  error  and  this  uncertainty  will  then  be 
compared  with  the  uncertainties  derived  from  the  experimental  volume 
determinations  shown  in  table  10.  Finally,  the  uncertainty  in  determining  V 
will  be  estimated. 

A-1-3-1  Mass  of  Water.  The  mass  of  water  required  to  fill  a cylinder  was 
determined  by  weighing  the  cylinder  empty  and  then  full  using  the  substitution 
weighing  method  on  a 50  kg  balance.  The  accuracy  of  this  mass  determination 
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is  influenced  by  the  variability  of  the  balance,  the  calibration  of  the 
weights,  and  the  correction  for  buoyancy  of  the  cylinder  and  the  weights.  The 
mass  of  water  was  computed  using  the  equation 

M = - W.  + , (17) 


where 


M = mass  of  water,  grams, 

= sum  of  face  values  of  the  weights  used  to  balance  the  filled 
cylinder , grams , 

W-  ■=  sum  of  face  values  of  the  weights  used  to  balance  the  empty  open 
cylinder,  grams, 

Cy  = sum  of  the  corrections  applied  to  the  face  values  of  the  weight, 
grams 

= buoyancy  correction,  grams,  and 

= correction  for  the  inequality  of  balance  arm  lengths,  grams. 


A-1-3-1-1  Balance  and  Weights.  The  variability  of  the  50  kg  balance  was 
assessed  from  the  reciprocal  sensitivity,  with  a given  weight  on  both  pans, 
and  from  the  reproducibility  of  the  rest  point  for  that  weight.  With  nominal 
loads  on  the  pans  of  30.5  kg  and  60.1  kg,  corresponding  to  the  initial  and 
final  weights  of  a cylinder,  the  largest  reciprocal  sensitivities  observed 
were  76  mg/division  and  112  mg/division,  and  the  greatest  uncertainties  in  the 
rest  point  reproducibilities  were  0.6  division  and  0.8  division,  respectively. 
Therefore,  the  maximum  uncertainty  arising  from  the  nonreproducibility  of  the 
balance  was  46  mg  in  W-  and  90  mg  in  W^,  so  that  the  standard  deviations  were 
15  mg  and  30  mg,  respectively. 


The  balance  arm  length  ratio  as  experimentally  determined  by 
transposition  weighing  differs  from  unity  by  2.5  ppm.  This  introduced  a 
systematic  error  in  the  differential  cylinder  weight  of  (W^  - W-)  x 2.5  x 10'^ 
gram  or  74  mg.  A corresponding  correction  was  applied,  and  the  residual 
uncertainty  was  probably  no  greater  than  1/2  ppm  or  15  mg  so  that  the  standard 
deviation  of  was  5 mg.  The  1/2  ppm  uncertainty  contributes  maximum  errors 
of  30  mg  and  15  mg  in  and  respectively. 

The  value  for  was  computed  from  the  calibration  corrections  supplied 
by  the  NBS  Length  and  Mass  Division  for  the  specific  set  of  Class  S weights 
[3]  used,  which,  for  weights  of  10  grams  and  above  were  determined  with  a 
maximum  uncertainty  of  3 ppm.  Since  - W-  was  29.6  kg,  it  follows  that  the 
maximum  uncertainty  in  was  89  mg  and  the  standard  deviation  of  was  30 
mg.  Furthermore,  the  3 ppm  uncertainty  contributes  maximum  uncertainties  in 
the  corrected  values  of  and  W.  alone  of  180  mg  and  90  mg,  respectively. 

A-1-3-1-2  Buoyancy  Correction.  It  can  be  shown  that  the  buoyancy  correction, 
C,^,  is  given  by 
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Cb  = (W, 


W,) 


pf 

Ph  0 
2 


w. 


(18) 


where 

= sum  of  corrected  values  of  the  weights  used  to  achieve  equilibrium 
with  the  filled  cylinder,  grams, 

W-  = sum  of  corrected  values  of  the  weights  used  to  achieve  equilibrium 
with  the  empty  open  cylinder,  grams, 

p-  = density  of  the  air  at  the  time  of  the  weighing  of  the  empty 
cylinder,  g/cm^, 

= density  of  the  open  and  empty  cylinder,  g/cm^,  and 

Py  = density  of  the  weights,  g/cm^. 

p^  = density  of  the  air  at  the  time  of  the  weighing  of  the  filled 

cylinder,  g/cm^,  and 

Pj^  Q = density  of  the  water  in  the  cylinder,  g/cm^. 

2 

The  magnitude  of  Cb  is  about  31  grams.  Consider  now  the  errors  involved  in 
computing  Cb* 

The  air  density  was  calculated  using  eq.(7).  The  estimate  of 
uncertainty  in  the  air  density  is  s(p.)  = s(p^)  = 0.42  x 10'^  g/cm^. 

The  estimate  of  the  uncertainty  in  the  water  density  is  s(pj^  q)  = 

45  X 10'^  g/cm^.  The  error  in  the  assumed  density  of  the  weights  ^ay  be  as 
large  as  ±0.1  g/cm^,  giving  rise  to  a systematic  error  in  Cb  of  ±0.06  gram. 

The  cylinders  and  the  attached  valves  were  fabricated  from  stainless 
steel  which  has  a nominal  density  of  7.8  g/cm^.  However,  attached  to  the 
cylinders  were  a few  small  brass  fittings  used  in  the  filling  operation  and 
coupled  to  the  valves  were  aluminum  air  actuators . The  uncertainty  in  the 
nominal  density  is  systematic  and  may  be  as  great  as  ±0.2  g/cm^  leading  to  a 
maximum  contribution  to  the  buoyancy  error  of  ±0.001  gram. 


The  maximum  uncertainties  for  W.  due  to  the  balance  variability,  the 
calibration  error,  and  the  inequality  of  arms,  were  46,  90  and  15  mg  for  a 
total  of  151  mg;  for  the  corresponding  contributions  were  90,  180  and  30  mg 
for  a total  of  300  mg.  The  standard  deviations  will  be  assumed  to  be  s(W. ) = 
50  mg  and  s(W^)  = 100  mg. 


The  following  nominal  values  were  utilized  in  calculating  the  numerical 
values  of  the  partial  derivatives.  W-  = 30.55  kg,  = 60.13  kg,  p-  = p^  = 

1.161  X 10'^  g/cm^,  Ph2o  = 0.997  g/cm^,  p^  = 8.4  g/cm^,  p^^  - 7.8  g/cia^ . 

10"^  g/cm^. 


Pf  - Pi  = 
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The  steps  in  the  computation  of  the  uncertainty  in  due  to  random 
effects  are  given  in  table  11.  The  estimate  for  the  standard  deviation  s(C|^) 
is  0.113  gram. 

A-1-3-1-3  Error  in  Mass  Determination.  The  random  uncertainty  in  the  mass 
determination  is  computed  from  the  equation 


s(M)  - [s2(wp  + s2(Wj)  + s2(C„)  + s2(C^)  + s2(C|,)]1/2 


(19) 


Note  that  in  this  case,  and  W-  are  the  uncorrected  values  [as  opposed  to 
the  case  in  eq.  (18)],  and  their  errors  depend  only  on  the  balance  variability 
and  arm  ratio  errors.  The  maximum  errors  from  these  sources  are  120  mg  in 
and  61  mg  in  W. , leading  to  the  respective  standard  deviations  of  40  mg  and  20 
mg. 

The  calculations  in  table  12  show  that  s(M)  = 0.125  g. 

A-1-3-2  Water  Density.  The  sources  of  error  in  the  density  of  the  water 
used  in  the  cylinder  volume  determination  are  the  same  as  those  discussed  in 
A-1-3-1.  The  same  value  for  the  standard  deviation  of  the  density  estimated 

AT  ^ 

there  will  be  used  here,  that  is,  45  x 10’°  g/cm-^. 

A- 1-3-3  Cylinder  Valve  Correction.  The  cylinder  volume  is  that  volume 
enclosed  by  the  shut  inlet  valve,  barometer  valve,  and  vacuum  gage  valve.  In 
filling  the  cylinder  with  water  for  calibration,  a procedure  was  followed 
whereby  water  was  unavoidably  or  inadvertently  trapped  in  one  or  more  of  the 
valves.  This  trapped  water,  which  filled  the  channel  in  the  ball  of  the  valve 
was,  of  course,  weighed.  A volumetric  correction,  based  on  the  geometry  of 
the  channel,  was  therefore  made  to  reduce  this  systematic  error.  This 
correction  was  computed  by  the  equation 


(20) 


where  D is  the  diameter  of  the  cylindrical  channel  and  is  the  length  of  the 
channel.  Since  D is  3/8 -in  and  L = 3/4- in  for  the  1/4- in  valve  (water  was 

C ^ 

not  trapped  in  the  large  valve),  C^  is  1.36  cm-^. 

The  maximum  uncertainty  AC^  is  estimated  from  the  equation 


5C, 


V 


ac 

V 


AC, 


V 


AD  + ^ AL, 


(21) 
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where 


aC^  nDL^ 
dD~  ^ 2 


(22) 


(23) 


AD  = maximum  uncertainty  in  D,  and 
AL^  = maximum  uncertainty  in  L^. 

At  the  time  this  work  was  done , normal  engineering  tolerances  permitted 
a maximum  error  of  ±1/64- inch  in  any  dimension.  In  the  absence  of  any  direct 
measurements  of  D or  L^,  it  will  be  assumed  that  the  maximum  uncertainty  in 
each  is  ±1/64- in  from  which  it  follows  that  AC„  for  a single  valve  is  ±0.14 
cm^  and  for  three  values,  ±0.42  cm-^.  This  error  is  systematic.  The  standard 
deviation  s(V^)  as  indicated  in  table  13,  is  1.34  cm^. 

A- 1-3-4  Correction  to  a Standard  Temperature.  The  following  correction  was 
applied  to  each  experimentally  determined  cylinder  volume  in  order  to  adjust 
it  to  the  value  it  would  have  at  the  standard  temperature  of  25‘’C: 


Cx  - - t^)  . (24) 

where 

= correction,  cm^, 

a = coefficient  of  volumetric  expansion  for  the  cylinder, 
cm^/  ( cm^  • * C ) , 

t^  = standard  temperature,  i.e.,  25*0, 

tg  = temperature  at  which  volume  V was  determined,  ®C,  and 
Vq  = volume  of  cylinder  at  temperature  t^,  cm^. 

The  cylinder  was  fabricated  from  type  316  stainless  steel. 

Interpretations  of  the  data  of  Furman  [20] , Beenaker  and  Swenson  [21] , Lucks 
and  Deem  [22] , and  the  Metals  Handbook  [23]  lead  to  values  of  the  coefficient 
of  linear  expansion  of  316  stainless  steel  of  15 . 2 , 15.6,  15.3,  and  15.8,  all 
in  units  of  10‘^  cm/(cm- °C) . If  all  these  values  are  given  equal  weight,  the 
mean  is  15.48  x 10'^  cm/cm/'’C.  Assuming  that  the  coefficient  of  volumetric 
expansion  is  three  times  the  coefficient  of  linear  expansion,  then  the  mean  a 
= 46.4  X 10'^  cm^/(cm^* °C)  with  a standard  deviation  of  the  mean  of  0.45  x 
10'^  cm^/(cm^-°C.  The  correction  per  °C,  that  is,  C^/(tQ-tg),  is  29700  cm^  x 

46.4  X 10'^  cm^/(cm^*  °C)  or  1.38  cmv°C  where  29700  cm^  is  the  nominal  cylinder 
volume . 
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The  following  nominal  values  of  the  parameters  were  used  in  the 
calculation: 

a = 46.4  X 10'^  cm^/(cm^- °C)  , 

= 29700  cm^,  and 
(to-t,)  = 2»C. 

The  computations  and  results  are  shown  in  table  14.  The  standard 
deviation  s(C^)  is  0.14  cm^. 

A-1-3-5  Random  Error  in  the  Measurement  of  Gas  Volume.  The  estimated 
standard  deviation  of  the  volume  adjusted  to  25°C,  s(V^),  is  given  by  [s^(V^) 

+ or  1.35  cm^.  The  estimated  standard  deviation  of  the  mean  volume 

s(V|^),  based  on  n repeated  determinations  of  V^,  is  s(V^)//n  which  reduces  to 
0.60  cm^  for  n = 5.  The  experimental  standard  deviations  of  the  mean  volume 
determinations  were  2.15  cm^  and  0.99  cm^  for  cylinder  Nos.  1 and  2 
respectively  (table  10) . It  is  apparent  that  the  estimate  of  the  uncertainty 
in  the  volume  determination  based  on  an  error  analysis  yields  a value  that  is 
of  the  same  order  of  magnitude  as  the  experimental  uncertainty. 

If  the  gas  volume  measuring  system,  during  a moisture  determination,  is 
maintained  at  some  temperature  other  than  25°C,  then  the  mean  volumes,  listed 
in  table  10,  must  be  adjusted  to  the  test  temperature  t by  adding 


a(t  - 25‘*C)  , (25) 

as  indicated  in  eq.  (15). 

It  follows  that  for  each  degree  difference  between  t and  t^,  the 
standard  deviation  of  is  0.02  cm^  and  therefore, 

s(V)  - [s2(V^)  + 0.0004(t-25')2]'/2  . (26) 

A-1-3-6  Systematic  Errors  in  the  Gas  Volume.  The  non-negligible  systematic 
errors  in  the  cylinder  volumes  and  thus  in  the  gas  volumes  arise  in  the 
buoyancy  correction  (A-1-3-1-2)  and  the  cylinder  valve  correction  (A-1-3-3). 
The  former  leads  to  a maximum  uncertainty  in  and  thus  in  the  water  mass  of 
±0.061  gram.  This  corresponds  to  a maximum  uncertainty  in  the  cylinder  volume 
of  ±0.06  cm^.  The  latter  leads  to  a maximum  uncertainty  in  the  cylinder 
volume  of  ±0.42  cm^.  Altogether  there  is  a maximum  possible  systematic  error 
in  the  cylinder  volume  of  ±0.48  cm^. 

A-1-4  DETERMINATION  OF  THE  GAS  DENSITY  IN  A CYLINDER  (see  [1]  pp. 418-423) 

The  gas  of  principal  interest  in  humidity  measurement  is  atmospheric 
air;  hence,  the  discussion  and  error  analysis  will  be  confined  to  this  gas 
although  similar  analyses  may  be  made  for  other  gases. 
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A-1-4-1  Standard  Density.  Air  is  a mixture  of  oxygen,  nitrogen,  argon  and 
carbon  dioxide.  The  percentage  composition  of  the  components  is  essentially 
constant.  Such  variables  as  geographic  location,  altitude,  and  presence  of 
sources  of  contamination  have  little  detectable  effect  on  the  percentage 
composition  [24].  The  percentages  by  volume  as  given  by  Harrison  [25]  are  as 
follows:  Nitrogen,  78.084;  Oxygen,  20.948;  Carbon  Dioxide,  0.0314  and  Argon, 
0.934.  These  values  are  based  generally  on  the  average  of  results  of  chemical 
analyses  made  at  various  times  and  places  [25]. 

The  value  of  p^,  the  density  of  air  at  273. 16K  and  760  mm  Hg,  is  computed 
from  the  equation 


M. 


M. 


V,  Z V 
a o 


(27) 


where  is  the  molecular  weight  of  air  in  g/mole , and  is  the  molar  volume 
in  cm^/mole  that  air  occupies  under  standard  conditions.  Since  air  is  not  an 
ideal  gas,  is  equal  to  V,  the  molar  volume  an  ideal  gas  occupies  under 
standard  conditions,  multiplied  by  Z^,  the  compressibility  factor  for  air  at 
273. 16K  and  760  Torr. 

The  molecular  weight  of  air  is  given  by 


Ma  = %N2(Mjj  ) + %02(Mq  ) + %C02(Mj.o  ) + %A(M^) 
2 2 2 


(28) 


where  the  M's  and  their  corresponding  subscripts  denote  the  molecular  weights 
of  the  four  components  listed  above. 


Therefore  is  given  by 

p^  = [ %N2(M^  ) + %02(Mq  ) + %C02(Mj.o  ) + %A(M^)  ]/VZ^  . (29) 

2 2 2 

Using  the  above  percentages  for  the  components  and  the  following  value 
for  their  molecular  weights,  N2  = 28.016  g/mole,  O2  = 32  g/mole,  CO2  = 44.011 
g/mole,  and  A = 39.944  g/mole,  and  using  V = 22414.6  cm^  [26],  and  Z^  = 
0.99941[4],  it  is  found  that  p^  = 1.29304  x 10'^  g/cm^. 

The  fact  that  there  are  uncertainties  in  p arising  from  uncertainties  in 
composition  will  be  taken  into  account  by  assuming  these  uncertainties 
contribute  to  the  standard  deviation  of  p^. 

According  to  information  available  at  the  National  Institute  of 
Standards  and  Technology  (formerly  National  Bureau  of  Standards)  [27] , the 
contribution  of  maximum  uncertainties  in  the  molecular  weights  of  the 
constituents  to  the  maximum  uncertainty  of  the  molecular  weight  of  air  is  4 or 
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5 parts  in  10^  from  nitrogen,  and  a maximum  of  1 part  in  10^  for  the  remaining 
components.  This  gives  a total  maximum  uncertainty  of  6 parts  per  10^. 
Applying  the  3 sigma  rule,  the  standard  deviation  of  the  molecular  weight  of 
air  is  2 parts  in  10^  . The  molecular  weight  of  air  is  28.966  g/mole  as 
computed  from  eq.  (28).  Therefore  its  standard  deviation,  arising  from 
uncertainties  in  the  molecular  weights  of  the  components,  is  6 x 10'^  g/mole. 

The  standard  deviation  of  V is  0.9  cm^  [26]. 

According  to  Hilsenrath  [4],  the  uncertainty  in  is  nominally  zero. 

From  the  computations  in  table  15,  the  standard  deviation  of  is  0.1  x 
10'^  g/cm^. 

The  density  of  air  at  any  other  temperature  and  pressure  is  given  by 


273.15  B 

0 

T ^ 7^  ^ z”  ’ 


(30) 


where  Z,  the  compressibility  factor,  is  a function  of  T and  B. 

Hilsenrath  [4]  has  shown  that  the  differences  between  compressibility 
factors  of  his  tables  and  the  tables  obtained  by  other  researchers , at  a 
nominal  pressure  of  1 atm,  is  zero  at  273. 16K  and  1 part  in  10^  at  323. 15K.  If 
this  is  interpolated  linearly,  the  difference  increases  2 parts  in  lOV^. 
Therefore  at  300K,  the  difference  is  5.4  parts  in  10^.  This  will  be  taken  as 
the  standard  deviation  of  the  compressibility  factor  for  the  present  purpose. 

The  compressibility  factor  at  one  atmosphere  and  300K  is  nominally 
unity,  so  its  standard  deviation  is  54  x 10'^. 

A- 1-4- 2 Pressure . The  accuracy  with  which  the  pressure  in  a cylinder  can  be 
measured  depends  on  the  barometer  errors,  the  residual  gas  pressure  at  the 
cessation  of  evacuation,  the  fluctuations  at  the  assumed  equilibrium  point, 
leakage  in  the  barometer  system  and  the  error  due  to  the  difference  between 
the  pressure  in  the  barometer  system  and  that  in  the  cylinder  at  the  moment 
the  barometer  valve  is  opened. 

A precision  micrometer  cistern  mercury  barometer  is  used  to  measure  the 
equilibrium  pressure  within  a cylinder.  This  instrument  was  calibrated  by  the 
NES  Pressure  and  Vacuum  Section  in  the  latter  part  of  1960.  After  applying 
corrections  for  the  zero  reading,  capillary  depression,  temperature,  and 
gravity,  it  is  estimated  that  the  standard  deviation  of  a pressure 
measurement,  arising  from  the  barometer  itself,  is  0.08  Torr. 

The  barometer  scale  was  independently  calibrated  with  gage  blocks,  with 
another  mercury  manometer,  and  with  a piston  gage.  At  10® C and  700  Torr,  the 
gage  block  calibration  indicates  a scale  correction  of  -0.005  Torr,  whereas 
both  of  the  other  methods  indicate  nominal  scale  corrections  of  +0.120  Torr. 
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Because  of  this  discrepancy,  a correction  of  +0.065  Torr  will  be  applied 
to  any  reading,  and  a residual  systematic  error  of  ±0.075  Torr  will  be 
assigned  to  any  reading. 

The  equilibrium  values  of  the  pressures  in  a given  cylinder  with 
repeated  evacuation  and  filling  have  varied  as  much  as  1.3  Torr,  although  a 
more  typical  value  would  be  0.5  Torr.  The  variations  arise  from  the  photocell 
circuit  response  time  and  the  variability  of  the  flow  rate  at  this  part  of  the 
cycle.  Because  of  the  magnitude  of  the  variations,  the  actual  pressure  must 
be  recorded  each  time  equilibrium  has  been  reached. 

Fluctuations  also  occur  after  the  cylinder  has  presumably  attained 
equilibrium.  Since  it  is  not  known  whether  these  fluctuations  arise  from  the 
nonreproducibility  of  the  barometer  or  are  real  fluctuations,  the  latter  will 
be  assumed,  in  which  case  this  contributes  another  uncertainty  to  the 
pressure.  The  largest  such  fluctuation  has  been  about  0.2  Torr,  although  the 
more  usual  figure  is  0.025  Torr.  The  standard  deviation  will  be  assumed  to  be 
represented  by  1/3  of  the  maximum,  that  is,  0.07  Torr. 

The  total  standard  deviation  of  a single  pressure  reading,  because  of 
errors  in  the  barometer  itself  and  the  fluctuations  which  occur  is  given  by 

s(p)  = [(0.08)2  + (0.07)2]1/2  . (31) 

or  s(p)  =0.10  Torr. 

At  least  three  barometer  readings  are  taken  after  the  cylinder  has 
reached  its  equilibrium  pressure,  and  these  are  averaged.  The  standard 
deviation  of  the  average  pressure  within  a cylinder  is  thus  given  by 


s(p) 

0.06  Torr.  (32) 


The  barometer  system  (i.e.,  the  cistern  and  connecting  tubing)  is  not 
necessarily  at  the  equilibrium  pressure  of  the  cylinders  before  a run  is 
started.  As  a result,  air  from  the  system  is  dumped  into  a cylinder,  or  vice 
versa  (depending  on  which  is  initially  greater) , resulting  in  a pressure 
reading  that  depends  both  on  the  true  cylinder  pressure  and  on  the  pressure  in 
the  barometer  system  before  connection  to  the  cylinder.  The  excess  air  is 
shunted  back  and  forth  between  cylinders  as  they  open  and  close,  so  that  the 
correction  need  only  be  applied  once,  upon  completion  of  the  run. 

The  correction  is  of  the  form 


VgAP  + A(P^^  - pf) 


(33) 
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where 


Vg  = barometer  system  volume  when  column  and  cistern  are  at  same 
pressure,  cm^, 

= final  cylinder  pressure  reading  on  barometer,  cm, 

P.  = initial  pressure  reading  on  barometer  before  run  begins,  cm, 

AP  = P^  - Pp  cm, 

A = area  of  barometer  column,  cm^,  and 

= volume  of  cylinder  to  which  the  barometer  is  finally  connected, 
cm^ . 

The  same  correction  applies  for  the  pressure  switch,  although  for  it  P^ 
and  P.  represent  final  and  initial  ambient  pressures.  The  sum  of  the  two 
corrections  is  about  0.1  Torr. 

The  residual  errors,  after  applying  the  correction  to  both  systems,  are 
estimated  to  be  one  part  in  40,000  if  the  pressure  reading  is  100  Torr  or 
greater,  that  is,  when  a cylinder  has  collected  1/7  of  its  volume,  or 
nominally  4.2  liters. 

The  leak  rate  of  the  barometer  system  is  such  that  the  mercury  column 
moves  0.006  mm/min  (essentially  equivalent  to  .006  Torr/min)  under  a vacuum  of 
less  than  100  millitorr. 

The  length  of  time  the  barometer  system  is  open  to  a cylinder  is  roughly 
10  minutes,  the  latter  half  of  which  is  under  equilibrium  conditions.  If  the 
leak  rate  was  0.006  Torr/min  at  the  equilibrium  pressure  of  695  Torr,  the 
column  would  be  expected  to  move  upward  0.03  Torr  during  the  reading  time. 

This  is  not  observed,  indicating  that  error  from  the  leak  rate  at  the 
equilibrium  pressure  is  negligible. 

Each  cylinder  was  repeatedly  evacuated  to  determine  its  residual 
pressure  just  prior  to  filling.  These  tests  showed  that  each  cylinder  reached 
50  millitorr  in  about  93  seconds. 

At  the  time  of  the  momentary  closing  and  reopening  of  the  vacuum  valve 
of  the  cylinder  being  evacuated,  i.e.,  120  seconds  after  the  start  of 
evacuation,  the  cylinder  pressure  was  about  37  millitorr.  The  closing  and 
reopening  of  the  valve  is  a deliberate  action  to  remove  gas  trapped  in  the 
valve  channel.  It  causes  the  cylinder  pressure  to  rise  again  to  approximately 
75  millitorr  after  which  it  returns  to  the  value  it  had  prior  to  the  closing 
in  approximately  90  seconds . 

At  the  maximum  flow  rate  of  2 fpm  of  air  into  one  cylinder,  the  other 
cylinder  was  evacuated  for  about  420  seconds,  at  which  time  the  pressure  in 
that  cylinder  was  reduced  to  about  23  millitorr. 

Both  cylinders  had  residual  pressures  of  less  than  20  millitorr  after 
evacuation  for  600  seconds,  and  after  2100  seconds  the  pressures  had  decreased 
to  17  millitorr. 
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The  residual  back  pressure  within  the  cylinder  when  it  is  opened  for 
filling  after  evacuation,  as  measured  with  the  system's  thermocouple  vacuum 
gage,  is  no  greater  than  20  millitorr.  This  is  a systematic  error  for  which  a 
correction  may  be  applied.  However,  the  magnitude  of  this  error,  about  1 part 
in  35,000,  is  so  small  that  it  may  be  neglected.  The  scatter  of  the 
observations  after  repeated  evacuations  is  about  10  millitorr,  which  is  only  1 
part  in  70,000  of  the  equilibrium  pressure  of  695  Torr  and  thus  negligible. 
Therefore,  the  only  random  factor  contributing  to  the  accuracy  of  the  cylinder 
pressure  that  is  significant  is  the  barometer  error  which  is  0.06  Torr. 

A- 1-4- 3 Temperature . The  uncertainties  in  the  temperature  measurement  of  the 
gas  in  the  cylinder  arise  from  the  error  in  the  thermopile  calibration,  the 
accuracy  of  the  potentiometer  used  for  measuring  the  thermopile  voltage,  and 
the  fluctuations  about  the  mean  equilibrium  temperature. 

Each  thermopile  was  calibrated  at  nominal  temperatures  of  20,  25,  30  and 
35 °C  against  a NBS  calibrated  platinum  resistance  thermometer  for  which  the 
estimated  maximum  error  is  0.001°C.  The  thermopile  has  an  output  of  about  160 
pV/°C;  the  potentiometer  used  to  measure  this  output  has  a maximum  uncertainty 
of  1 pV;  therefore,  the  calibration  data  had  an  equivalent  maximum  uncertainty 
of  0.006®C.  In  order  to  permit  interpolation  between  thermopile  calibration 
points,  the  data  were  fitted  to  an  empirical  equation  of  the  form 

E = at  + bt^  + ct^  + dt^  . (34) 

This  equation  was  used  to  compute  the  EMF  outputs  for  temperatures  at  0.5‘’C 
intervals,  and  straight  lines  were  then  drawn  between  successive  points.  These 
lines  are  used  to  convert  EMF  readings  into  temperature.  A comparison  of 
interpolated  values,  as  derived  from  the  curves  and  the  experimental  data 
disclosed  that  the  maximum  difference  was  0.011°C.  Subsequent  to  the 
calibration,  each  measurement  with  the  thermopile  involves  the  use  of  the 
potentiometer  so  that  the  latter  again  contributes  an  uncertainty  of  0.006°C. 

The  equilibrium  temperature  fluctuated  about  0.002°C.  Since  it  is  not 
known  whether  this  is  because  of  the  nonreproducibility  of  the  thermopiles  or 
whether  it  is  real,  an  error  of  this  magnitude  will  be  assigned  to  the 
temperature . 

The  sum  of  the  uncertainties  in  the  measurement  of  cylinder  temperature 
is  thus  0.026®C.  It  will  be  assumed,  therefore,  that  the  standard  deviation 
of  the  gas  temperature  in  the  cylinder  at  equilibrium  is  0.009  deg  C,  that  is, 
1/3  of  the  maximum  uncertainty. 

A-1-4-4  Random  Error  in  the  Determination  of  Gas  Density.  Using  the 
appropriate  partial  derivatives  of  eq.  (29),  along  with  the  required  standard 
deviations , leads  to  the  terms  contributing  to  the  uncertainty  in  the  gas 
density.  The  nominal  values  used  for  the  independent  parameters  are  Pq  = 1.3 
X 10'^  g/cm^,  Z /Z  = 1,  T = 298K,  and  B = 700  Torr.  As  indicated  in  table  16, 
s(p)  = 0.14  X 10'°  g/cm^.  These  values  vary  slightly  with  temperature,  but 
the  variations  are  insignificant  for  the  present  purpose. 

A- 1-4 -5  Systematic  Error  in  the  Determination  of  Gas  Density.  The  only 
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non-negligible  systematic  error  occurring  in  this  determination  is  that 
arising  from  the  uncertainty  of  the  scale  correction  to  the  barometer.  As 
mentioned  in  A-1-4-2,  the  magnitude  of  this  error  is  ±0.075  Torr. 

Use  of  the  density  equation  [eq.  (7)]  shows  that  the  magnitude  of  the 
corresponding  systematic  error  in  the  air  density  is  ±0.12  x 10’^  g/cm^. 

A- 2 ANALYSIS  OF  COMPLETE  SYSTEM 

A-2-1  ACCURACY  IN  THE  DETERMINATION  OF  MIXING  RATIO  (see  [1]  pp. 423-428) 

A-2-1-1  Random  Error  in  the  Mixing  Ratio.  The  mixing  ratio  is  defined  by  eq. 
(1)  while  the  standard  deviation  of  the  mixing  ratio  is  given  by  eq.  (2).  The 
latter  equation  is  repeated  here. 


dr 

dr 

dr 

s2(m) 

s2(p))’/2  , 

(35) 

where 


dr  1 
dM  ■ Vp 


(36) 


dr 

dV 


and 
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^ = _M 
dp  Vp' 


(38) 


It  will  be  convenient  to  use  relative  errors  rather  than  absolute  errors 
at  this  time.  If  both  sides  of  eq.  (35)  are  divided  by  the  mixing  ratio  r, 
then 


s(r) 

r 


s^(M)  s^(V)  s^(p) 


1/2 


(39) 
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The  nominal  value  of  M is  dependent  on  the  duration  of  a test  run,  the 
flow  rate  of  the  test  gas,  and  the  moisture  content  of  the  test  gas,  whereas 
s(M)  is  essentially  independent  of  these  parameters.  If  the  test  gas  has  a 
high  moisture  content,  then  at  the  design  flow  rates  of  this  apparatus,  in  a 
relatively  short  time  enough  moisture  can  be  collected  in  the  drying  train  so 
that  the  relative  errors  s(M)/M  is  reasonably  small.  However,  as  the  moisture 
content  of  the  test  gas  decreases,  to  collect  the  same  mass  of  water  the 
duration  of  a run  must  be  increased.  At  a given  flow  rate  a low  enough 
moisture  content  will  eventually  be  reached  for  which  the  duration  of  a run 
will  be  inordinately  long  to  maintain  the  same  relative  error.  Thus  for 
moisture  contents  below  this  value,  the  duration  of  the  run  will  be  the  factor 
limiting  the  accuracy. 

It  has  been  estimated  that  M can  be  determined  with  a standard  deviation 
of  13  X 10'^  grams  provided  all  the  moisture  is  removed  by  the  first  two 
U- tubes  in  the  main  drying  train.  Assuming  that  0.6  gram  of  water  are 
collected,  the  relative  error  is  thus  13  x 10’^  /0.6  or  2.2  parts  in  10'^. 

The  volume  of  the  test  gas  will  be  the  sum  of  the  nuimber  of  times  each 
cylinder  is  filled  multiplied  by  its  calibrated  volume.  The  error  is  the  same 
for  each  cylinder,  so  that  the  total  error  is  the  number  of  fillings  times  the 
error  for  a single  filling.  Hence  the  relative  error  in  the  total  dry  air 
volume  is 


V 


(40) 


irrespective  of  the  number  of  times  a cylinder  is  filled. 

As  indicated  in  A-1-3-5,  the  error  is  a function  of  temperature.  At 
25° C,  the  experimental  standard  deviations  of  the  mean  volumes  were  2.15  and 
0.99  cm^  for  cylinders  1 and  2,  respectively,  whereas  the  computed  value  was 
0.60  cm^.  The  value  2.15  cm^,  being  the  largest  of  the  three,  will  be  used  as 
the  estimate  of  the  standard  deviation  for  the  cylinder  volume. 

Consider  now  eq.  (25)  where  the  experimental  standard  deviation  of  the 
mean  value  is  used  instead  of  the  calculated  standard  deviation  of  a single 
volume  determination.  The  equation  becomes 

s(V)  = [(2.15)2  + o.0004(t  - 25°C)2]1/2  (41) 

The  maximum  value  of  t-25°C  is  10°C  for  the  instrument.  Then  s(V)  can 
be  shown  to  be  2.16  cm^.  The  relative  error  is  2.16/(29.7  x 10^)  or  0.73 
parts  in  10^.  Since  the  volume  error  is  almost  independent  of  the  bath 
temperature,  the  value  0.73  parts  in  10^  will  be  used  for  further  calculation. 
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The  nominal  value  of  p is  1.1  x 10’^  g/cm^.  The  estimated  s(p)  is  0.14 
X 10'^  g/cm^.  The  relative  uncertainty  in  p , i.e.,  s(p)/p,  is  therefore  1.3 
parts  in  10^. 

The  relative  standard  deviation  of  the  mixing  ratio  is  given  by 
s(r)  = [(1.3)^  + (0.73)^  + (2.2)^]^^^  , in  units  per  10^  where  1.3  = 

s(p);  0.73  = s(V);  and  2.2  = s(m).  The  latter  is  based  on  collecting  0.6  gram 
of  water.  Thus  s(r)  = 2.7  parts  per  10^. 

A- 2 -1-2  Systematic  Errors  in  the  Mixine  Ratio.  In  addition  to  the  systematic 
errors  which  arise  in  the  measurements  of  water  vapor  mass,  cylinder  volume, 
and  dry  air  density,  there  are  three  other  known  non-negligible  systematic 
error  sources  associated  with  the  operation  of  the  gravimetric  hygrometer. 
These  are  the  incompleteness  of  water  vapor  absorption  by  the  U-tubes,  an 
irregular  sampling  flow  rate  coupled  with  a possibility  of  variation  in  the 
test  gas  moisture  content,  and  leakage  of  room  air  into  the  drying  train  and 
cylinders . 

Only  the  extremes  of  these  errors  can  be  estimated,  although  the  actual 
error  in  any  test  could  lie  anywhere  between  zero  and  the  extremes. 
Accordingly,  no  attempt  is  made  to  correct  for  these  errors. 

The  error  arising  from  two  systematic  maximum  errors  of  the  same  sign  is 
given  by  the  algebraic  sum  of  the  two.  On  the  other  hand,  systematic  maximum 
errors  of  opposite  signs  may  not  be  added  algebraically.  Rather,  systematic 
errors  of  opposite  signs  give  rise  to  an  error  band.  The  upper  and  lower 
error  limits  are  not  necessarily  the  same. 

A-2-1-2-1  Systematic  Errors  in  Mass  of  Water  Vapor.  Cylinder  Volume,  and  Dry 
Gas  Density.  It  was  shown  in  A-1-2-9,  A-1-3-6,  and  A-1-4-5  that  the  maximum 
possible  systematic  errors  in  the  mass  of  water  vapor,  cylinder  volume  and  gas 
density  are  respectively  ±0.048  mg,  ±0.48  cm^,  and  ±0.12  x 10'^  g/cm^. 

Since  the  nominal  values  of  these  parameters  are  m = 0.6  gram,  V = 29.7 
liter,  and  p = 1.06  x 10'^  g/cm^,  the  relative  systematic  errors,  in  units  of 
parts  per  10^,  are  ±0.80,  ±0.16  and  ±1.13  (for  m,  V and  p) . 

Therefore,  these  sources  contribute  a maximum  systematic  error  in  the 
mixing  ratio  of  ±2.09  parts  per  10^. 

A-2-1-2-2  Incompleteness  of  Absorption  bv  the  U-tubes.  Although  the  data 
analysis  of  the  water  pick-up  experiments  (A- 1-2 -7)  revealed  no  systematic 
tendencies,  it  was  suggested  that  a maximum  of  0.5  x 10'^  mg  of  water  vapor 
may  go  undetected  for  each  liter  of  dry  air  collected. 

The  maximum  duration  of  the  water  pick-up  runs  was  nominally  8 hours  at 
a flow  of  about  1.5  Cpm.  Thus  720  liters  of  air  were  drawn  through  the 
absorption  tubes  and  the  mass  of  water  vapor  that  may  have  escaped  absorption 
was  at  most  0.5  x 10'^  x 720  or  0.04  mg.  If  in  any  run  the  amount  was  less, 
it  could  easily  have  gone  undetected. 
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Therefore,  the  experiment  described  in  A- 1-2-10  is  not  a conclusive 
basis  for  assuming  the  nonexistence  of  this  type  of  systematic  error, 
especially  for  tests  which  are  more  than  8 hours  long,  and  such  an  error  will 
be  assumed  to  be  present. 

The  nominal  air  density  encountered  during  tests  with  the  gravimetric 
hygrometer  is  1.06  g/C . Since  the  maximum  amount  of  water  that  may  be  in  the 
effluent  from  the  absorption  tubes  is  0.5  x 10'^  nig/{  of  air,  the  measured 
mixing  ratio  can  be  lower  than  the  actual  mixing  ratio  by  at  most  0.5  x 10'^ 
mg/1.06  gram  or  0.47  x 10'^  mg/g,  regardless  of  the  magnitude  of  the  mixing 
ratio . 


One  of  the  criteria  for  the  operation  of  the  gravimetric  hygrometer  is 
that  at  least  0.6  gram  of  water  vapor  shall  be  collected  in  any  test.  Thus  a 
test  made  at  a low  moisture  content  will  require  a larger  volume  of  test  gas 
to  be  drawn  through  the  U- tubes  and  therefore  more  cylinders  to  be  filled  than 
at  a high  moisture  content.  Assuming  that  the  mass  of  air  in  a filled 
cylinder  is  31.6  grams  and  that  0.6  gram  of  water  vapor  is  collected  per  run, 
a nominal  mixing  ratio  may  be  computed  for  a run  in  which  a given  number  of 
cylinders  are  filled. 

Table  17  lists  examples  of  the  number  of  cylinder  fillings  at  which  a 
run  may  be  terminated,  nominal  mixing  ratios  in  units  of  milligrams  of  water 
vapor  per  gram  of  associated  dry  air,  and  the  relative  error,  E,  of  the  mixing 
ratio  in  parts  per  10’^,  based  on  the  relation 


0.47  X 10'^ 
r 


E 

1^ 


(42) 


where  r is  the  mixing  ratio  in  mg/g. 

The  unabsorbed  moisture  passes  into  the  cylinders  where  it  is  measured 
as  dry  air.  Therefore  the  apparent  measured  amount  of  dry  gas  is  too  high  by 
the  amount  of  vapor  present,  and  this  again  causes  the  measured  mixing  ratio 
to  be  smaller  than  it  should  be. 

Since  the  molecular  weight  ratio  of  water  to  air  is  0.62197,  each  unit 
mass  of  water  vapor  displaces  1/0.62197  or  1.608  units  of  air  mass.  Thus  the 
0.5  X 10'^  gram  of  water  vapor  per  liter  of  dry  air  that  flows  into  the 
cylinder  causes  an  overstatement  of  the  amount  of  dry  air  present  by 
0.8  X 10'^  gram  for  each  liter  present. 

The  relative  error  X in  the  mass  of  dry  gas  is  found  from  the  relation 


0.8  X 10~^  ^ 

p 10^ 


(43) 


26 


where  p is  the  nominal  air  density  which  for  this  calculation  may  be 
considered  equal  to  1 g/H  . The  relative  error  is  thus  0.8  x 10’^  parts  in 
10^,  which  is  completely  negligible. 

A-2-1-2-3  Flow  Effect  Error.  At  average  sampling  flow  rates  between  1 {pm 
and  the  maximum  of  2 fpm,  the  rate  is  faster  when  a cylinder  begins  to  fill 
than  when  it  is  almost  full.  This  cyclic  sampling  affects  neither  the  flow 
through,  nor  the  output  of,  the  humidity  generator,  which  serves  as  the  source 
of  test  gas.  However,  if  the  output  of  moisture  content  from  the  generator 
were  to  be  cyclic  during  brief  time  periods,  then  the  possibility  arises  that 
the  sampling  rate  and  generator  output  could  be  synchronous  for  such  periods. 
Unless  the  total  length  of  the  run  were  long  compared  to  the  synchronous 
period,  a significant  systematic  error  in  the  average  value  of  the  mixing 
ratio  as  determined  by  the  gravimetric  hygrometer  could  result. 

For  example,  if  the  hygrometer  were  to  sample  low  moisture  content  air 
at  2 Cpm  for  10  minutes,  then  high  moisture  content  air  at  1 Cpm  for  10 
minutes,  it  would  sample  20  liters  of  low  moisture  content  and  only  10  liters 
of  high  moisture  content,  biasing  the  result  in  favor  of  the  former. 

Because  the  probability  of  synchronization  is  small,  the  averaging 
nature  of  the  operation  is  considered  to  make  this  flow  effect  error 
negligible  for  runs  in  which  the  number  of  cylinders  filled  is  greater  than 
three . 


For  runs  in  which  the  number  of  cylinders  filled  is  less  than  three,  the 
effect  may  be  eliminated  by  maintaining  the  flow  at  1 {pm  or  less.  At  these 
lower  rates,  the  flow  can  be  maintained  with  adequate  constancy.  Since  this  is 
the  manner  in  which  short  runs  are  performed,  the  flow  effect  is  negligible 
over  the  entire  operational  range  of  the  instrument. 

Subsequent  error  analysis  is  based  on  a flow  rate  of  1 {pm,  so  no 
additional  error  is  introduced  by  the  flow  restriction  for  runs  in  which  the 
number  of  cylinders  filled  is  less  than  three. 

A-2- 1-2-4  Leakage  Error.  A leakage  check  of  both  cylinders  indicates  a 
completely  negligible  leakage  rate  at  pressure  near  20  millitorr. 

When  a check  is  made  on  a system  consisting  of  one  cylinder  and  the 
pressure  switch,  the  leakage  rate  decreases  monotonically  from  0.03  cm^/min  in 
a pressure  range  of  20  to  50  millitorr  to  0.02  cm^/min  in  a pressure  range  of 
100  to  120  millitorr. 

At  a sample  gas  flow  rate  of  one  {pm,  a constant  leakage  rate  of  0.02 
cm^/min  would  contribute  an  extra  0.02  cm^  of  room  air  to  every  liter  of  dry 
test  gas,  an  error  of  two  parts  in  10^  in  the  volume  of  the  dried  test  gas. 
However,  since  the  leakage  rate  probably  decreases  below  0.02  cm/min  as  the 
pressure  increases,  the  error  it  contributes  to  the  air  volume  in  the 
cylinders  is  considered  to  be  negligible. 
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The  main  drying  train  was  subjected  to  separate  leakage  tests.  Air 
leaking  into  the  train  enters  from  the  room,  carrying  water  vapor  that  is 
absorbed  in  the  U- tubes.  The  tests  indicate  approximate  leakage  of  0.003 
cmVniin  at  pressures  in  the  millitorr  range,  and  0.002  cm^/min  at  a nominal 
pressure  of  340  Torr.  This  leakage  apparently  does  not  decrease  appreciably 
with  increasing  pressure,  and  it  will  be  assumed  that  the  leakage  rate  of 
0.002  cm^/min  holds  over  the  entire  range  of  pressures.  Furthermore,  it  will 
be  assumed  that  the  entire  leakage  occurs  into  the  first  two  tubes  of  the  main 
drying  train,  where  the  water  vapor  in  the  incoming  room  air  is  picked  up  and 
later  weighed  as  part  of  the  total  water  mass.  The  dried  air  from  the  leakage 
then  passes  on  into  the  cylinders,  where  it  is  measured  as  part  of  the  total 
dry  air  mass. 

The  mixing  ratio  is  given  by  eq.  (1).  However,  since  there  is  leakage, 
the  mixing  ratio  actually  measured  is  given  by 

r = M + p„3VV(Vp  + V'p')  , (44) 

where 

M,V,p  = the  quantities  defined  for  eq.  (1), 

pyg  = mass  of  water  vapor  per  unit  volume  of  dry  room  air,  g/H  , 

p'  = density  of  dry  room  air,  g/H  , and 

V'  = volume  of  dry  room  air  entering  the  system,  C . 

The  mass  of  dry  air  collected  per  cylinder,  Vp,  is  about  31.6  grams  and  the 
mass  of  water  collected  per  run,  M,  is  about  0.6  gram. 

If  a flow  rate  of  one  {pm  is  assumed,  then  it  takes  30  minutes  to  fill  a 
cylinder,  and  the  total  volume  of  room  air  leaking  into  the  drying  train 
during  the  filling  of  each  cylinder  is  6 x 10'^  liters.  Assuming  the  nominal 
room  air  density  to  be  1.17  g/H  , the  mass  of  dry  air,  per  cylinder,  from  the 
leakage  into  the  absorption  tubes  is  1.17  x 6 x 10'^  gram  or  7.02  x 10'^  gram. 

To  assess  the  quantity  PygV'  , the  water  vapor  mass  from  the  room  air 
that  is  collected  in  the  U-tubes,  it  is  necessary  to  consider  conditions  of 
both  maximum  and  minimum  water  vapor  density  in  the  ambient  air  in  the 
vicinity  of  the  drying  train. 

The  lowest  temperature  in  the  area  in  which  the  hygrometer  is  used  may 
be  about  19®C.  A room  relative  humidity  of  10  percent  is  assumed  possible  at 
this  temperature.  At  this  temperature,  the  saturation  concentration  of  water 
vapor  in  air  is  16.31  x 10’^  gram  of  water  vapor  per  liter  of  air,  so  that  at 
a relative  humidity  of  10  percent,  the  concentration  is  1.631  x 10’^  gram  of 
water  vapor  per  liter  of  air.  Since  6 x 10'^  liters  of  room  air  leak  into  the 
system  per  cylinder,  1.631  x 6 x 10'®  or  9.8  x 10'®  gram  of  water  vapor  are 
collected  from  the  leakage  for  each  cylinder  filling. 

The  highest  ambient  temperature  is  28° C,  with  a possible  relative 
humidity  of  60  percent.  The  saturation  concentration  is  27.24  x 10’^  gram  of 
water  vapor  per  liter  of  air,  so  that  at  a relative  humidity  of  60  percent, 
the  concentration  is  16.344  x 10'^  gram  of  water  vapor  per  liter  of  air.  Thus 
6 X 10'^  X 16.344  X 10'^  or  98.1  x 10'®  gram  of  water  vapor  is  collected  from 
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the  leakage,  for  each  cylinder. 

In  table  18  are  tabulated  the  masses  of  air  and  water  gained  form  both 
the  test  gas  and  the  leakage  gas  under  the  conditions  of  maximum  and  minimum 
observed  water  vapor  density  in  the  ambient  air,  for  tests  involving  different 
numbers  of  cylinder  fillings. 

Table  19  is  a continuation  of  table  18  in  which  the  "true"  mixing  ratio 
(based  on  the  assumptions  of  collecting  0.6  gram  water  vapor  per  run  and  that 
the  mass  of  associated  dry  air  per  cylinder  is  31.6  grams),  and  the  measured 
mixing  ratio  (assuming  leakage  in  addition  to  the  above)  are  given,  along  with 
the  relative  systematic  error. 

The  mixing  ratios  have  been  computed  to  as  many  figures  as  was  necessary 
to  get  an  error  indication.  The  ratios  are  listed  in  terms  of  mg/g. 

When  the  room  air  mixing  ratio  is  less  than  the  test  gas  mixing  ratio,  the 
error  sign  is  negative,  and  for  the  inverse  situation,  the  error  sign  is 
positive . 

A- 3 OVERALL  ERROR  IN  MIXING  RATIO 

The  discussion  of  systematic  errors  in  the  mixing  ratio  was  based  on 
estimates  of  the  maximum  contributions  to  the  systematic  uncertainty  in  a 
mixing  ratio  determination  expected  from  various  parameters. 

It  was  shown  in  A- 2 -1-1  that  random  errors  gave  rise  to  a relative 
standard  deviation  of  the  mixing  ratio  of  2.7  parts  in  10^.  The  maximum 
contribution  of  the  random  errors  to  the  uncertainty  in  a mixing  ratio 
determination  is  considered  to  be  three  standard  deviations  or  ±8.1  parts  in 
10^. 


Examples  of  runs  in  which  a given  number  of  cylinders  are  filled  are 
listed  in  table  20,  along  with  the  corresponding  nominal  mixing  ratios,  and 
the  maximum  relative  errors  and  their  sources. 

For  any  run,  the  maximum  negative  and  positive  limits  of  error  may  be 
determined  by  summing  separately  the  negative  errors  and  the  positive  errors. 
Consider  a run  in  which  only  0.7  of  a cylinder  is  filled,  as  an  example.  The 
negative  errors  listed  in  table  20,  in  units  of  parts  in  10^  were  -8.1,  -2.09, 
-0.017  and  -0.02,  the  sum  of  which  is  nominally  -10.2.  This  defines  the  lower 
error  limit  in  this  case.  The  positive  errors  are  8.1  and  2.09,  the  sum  of 
which  is  nominally  +10.2.  This  defines  the  upper  error  limit. 

Therefore,  for  a run  in  which  only  0.7  of  a cylinder  is  filled,  the 
error  in  mixing  ratio  will  lie  in  the  band  from  -10.2  to  +10.2  parts  in  10^. 
Similar  calculations  may  be  made  for  runs  involving  any  number  of  cylinder 
fillings.  A synopsis  of  such  calculations  is  presented  in  table  21. 

A-4  CONCLUSIONS  AND  DISCUSSIONS 


The  National  Institute  of  Standards  and  Technology  (formerly  National 
Bureau  of  Standards)  has  developed  a hygrometer  for  the  measurement  of  mixing 
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ratios  lying  in  the  range  from  27  to  0.19  mg/g. 

The  highest  measurable  value  is  determined  by  consideration  of  operator 
comfort.  An  ambient  room  temperature  of  slightly  above  30 °C  is  tolerable  for 
short  periods,  and  to  prevent  condensation  in  the  line,  the  room  temperature 
must  be  above  the  dew-point  temperature.  This  limits  the  upper  dew  point  of 
the  test  gas  to  30° C,  corresponding  to  a mixing  ratio  of  nominally  27  mg/g, 
although  the  hygrometer  itself  is  capable  of  measuring  a higher  value. 

The  lower  limit,  0.19  mg/g,  is  determined  primarily  by  error 
considerations,  since  it  is  seen  that  systematic  errors  from  leakage  and 
incomplete  absorption  become  important  in  this  range. 

The  estimated  maximum  error  magnitude  encountered  over  the  operational 
range  of  the  Instrument  is  12.7  parts  in  10^,  or  0.13  percent  of  the  measured 
value . 


B.  NBS  TWO-PRESSURE  HUMIDITY  GENERATOR 

A precision  humidity  generator  [28]  is  used  as  the  principal  facility 
for  calibrating  transfer  and  secondary  standards,  and  for  testing  and 
evaluating  hygrometers  and  sensors . The  generator  operates  on  what  is  now 
known  as  the  two-pressure  principle.  A stream  of  air  at  an  elevated  pressure 
is  saturated  with  respect  to  the  liquid  or  solid  phase  of  water  and  then 
expanded  to  a lower  pressure.  Measurements  of  the  pressure  and  temperature  of 
the  saturated  air  stream,  and  in  the  test  chamber  after  expansion,  yield  the 
data  necessary  to  compute  the  water  vapor  content  of  the  air  stream. 

The  calibration  of  a hygrometer  with  the  two-pressure  humidity  generator 
can  be  made  in  various  units  which  relate  to  the  quantity  of  water  vapor  in  a 
moist  air.  Among  the  most  common  units  are  mixing  ratio,  dew-point 
temperature,  relative  humidity,  and  volume  ratio.  These  expressions  are 
defined  in  terms  of  real  gas  behavior  and  account  for  the  fact  that  the 
saturation  pressure  of  pure  water  in  the  presence  of  an  inert  gas  differs  from 
that  of  pure  water  alone . 

B-1  COMPUTATIONS  OF  HUMIDITY 

The  saturation  mixing  ratio,  r^,  of  the  moist  gas  emerging  from  the 
generator  is 


V(I’s'Ts)e„(Ts) 


Mg[Ps-f(Ps-Ts)e«(T3)] 


(45) 
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where 


My  = molecular  weight  of  water  vapor, 

Mg  = molecular  weight  of  the  carrier  gas, 

ey(Tg)  = saturation  vapor  pressure  over  a plane  surface  of  the  pure 
phase  of  liquid  or  solid  water  at  the  saturator  temperature, 
T 

■‘‘s  ’ 

Pg  ^ saturator  pressure,  and 

f = enhancement  factor  and  is  defined  below. 

The  enhancement  factor,  f,  at  the  saturator  pressure,  P^,  and 
temperature,  T^,  is  expressed  by 


“ e (T  ) 
w'^  S' 


(46) 


where  Xg,  = the  mole  fractions  of  carrier  gas  and  water  vapor  in  the 
saturated  mixture,  respectively. 

The  definition  of  the  relative  humidity  RH  in  percent  in  the  test 
chamber  of  the  generator  is  expressed  by 

RH  = (Xv/Xy)pj.,^j.  X 100  , (47) 

where  x^  = the  mole  fraction  of  water  vapor  in  a given  sample  of  moist  air 
characterized  by  pressure,  P^,  and  temperature,  T^,  and  x^  = the  mole  fraction 
of  water  vapor  in  the  saturated  mixture  at  the  same  values  of  pressure  and 
temperature . 

Substituting  appropriate  expressions  for  the  mole  fractions  yields 


RH 


f(Pe.T,) 


X — X 100  , 


(48) 


where  f(P^,T^)  = enhancement  factor  at  test  chamber  pressure,  P^,  and 

temperature  T^, 

ey(Tc)  = saturation  vapor  pressure  over  a plane  surface  of  the  pure 
phase  of  liquid  or  solid  water  at  the  test  chamber 
temperature , T^ . 

The  " thermo d3mamic  dew-point  (or  frost-point)  temperature”  T^  of  a moist 
gas  at  absolute  total  pressure,  P,  is  defined  as  that  temperature  at  which  the 
moist  gas  is  saturated  with  respect  to  a plane  surface  of  pure  liquid  (or 
solid)  water.  The  dew  point,  T^,  of  the  moist  gas  of  the  two-pressure 
generator  is  obtained  by  the  iterative  solution  of 
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f(Pc,T^)%(Td)  = f(Ps«Ts)ew(Tg)Pc/Ps 


(49) 


where  e^  and  f are  values  obtained  from  suitable  tables  and  is  the  absolute 

w c 

pressure  in  any  space  or  volume  that  is  filled  with  the  moist  gas,  e.g.,  the 
mirror  chamber  of  a dew-point  hygrometer. 

The  volume  ratio,  V,  of  the  moist  gas  of  the  two -pressure  generator  is 


f(P,,T,)e„(T3) 


(50) 


Hence,  with  the  establishment  of  constant  temperature  and  pressures  in 
the  saturator  and  test  chamber  the  various  units  of  humidity  can  be  calculated 
for  the  moist  gas  produced  by  the  two -pressure  humidity  generator.  The 
formulations  of  Wexler  [29,30]  are  used  for  obtaining  the  saturation  vapor 
pressure  of  water  and  ice.  When  air  is  used  as  the  carrier  gas,  the  values 
for  the  enhancement  factor,  f,  for  air,  given  by  Hyland  [31]  over  the 
temperature  range  -50°  to  90° C and  pressures  from  0.25  x 10^  to  10^  Pa  are 
used  in  the  above  equations.  Greenspan  [32]  has  obtained  a simplified 
equation  for  f which  can  be  easily  programmed  for  a computer  or  can  be 
calculated  with  the  aid  of  a programmable  pocket  calculator. 

Since  both  dew/frost  point  and  relative  humidity  can  be  expressed  with 
respect  to  either  ice  or  water,  appropriate  functions  for  the  enhancement 
factor  and  the  saturation  vapor  pressure  are  required  in  eqs.  (48)  and  (49) 
for  computations . 

Intercomparison  tests  were  made  with  the  NBS  standard  gravimetric 
hygrometer  [1]  over  a portion  of  the  generator's  operating  range.  The 
estimated  maximum  uncertainty  (three  standard  deviations)  is  0.2  percent  RH 
for  temperatures  between  0°  to  80° C which,  in  units  of  dew  point,  corresponds 
to  an  estimated  maximum  uncertainty  of  0.04°C  for  dew  points  from  -35  to  80° C. 

B-2  ANALYSIS  OF  COMPONENT  SYSTEMS 


As  indicated  in  eqs.  45  through  50,  the  calculations  of  the  various 
units  of  humidity  require  the  measurement  of  the  temperatures  and  pressures  of 
the  final  saturator  and  the  test  chamber. 

To  facilitate  the  automation  of  data  aquisition,  the  resistances  of  the 
calibrated  four- lead  standard  platinum  resistance  thermometer  and  the 
calorimetric  type  platinum  resistance  thermometer  are  measured  with  a bridge 
based  on  a design  by  Cutkosky  [33].  The  bridge  utilizes  an  inductive  ratio 
divider  and  requires  only  one  adjustment  for  balancing.  A built-in 
phase- sensitive  null-detector  easily  resolves  1 pQ  in  25  Q.  Small  deviations 
from  balance  are  recorded  continuously  with  an  analog  recorder  and/or  the  BCD 
output  of  a digital  voltmeter  is  recorded  on  a line  printer  at  any  preselected 
time  interval.  The  platinum  resistance  thermometers  were  calibrated  at  NBS  on 
the  International  Practical  Temperature  Scale  of  1968  and  subsequently  checked 
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from  time  to  time  at  the  triple  point  of  water.  It  is  estimated  that  the 
uncertainty  in  the  temperature  measurement  is  an  order  of  magnitude  more 
accurate  than  the  required  10  millidegree  Celsius  in  order  to  maintain  an 
accuracy  of  0.2%  for  most  of  the  mixing  ratio  range. 

The  pressures  are  measured  with  calibrated  fused  quartz  Bourdon  tube 
pressure  gages  equipped  with  BCD  outputs.  The  ranges  of  the  pressure  gages 
used  in  the  generator  are  0 to  0.21  MPa  (30.5  psia)  for  the  test  chamber  and  0 
to  0.69  MPa  (100  psia)  or  0 to  3.45  MPa(500  psig)  for  the  saturator.  These 
gages  are  periodically  calibrated  with  a dead  weight  piston  gage.  The 
accuracy  of  the  pressure  measurement  is  estimated  at  70  Pa.  This  has  been 
taken  into  account  in  the  overall  random  and  systematic  errors. 

B-3  ANALYSIS  OF  COMPLETE  SYSTEM 


Two  independent  approaches  were  used  to  evaluate  the  performance  of  this 
generator.  First,  an  intercomparison  was  made  between  the  generator  and  the 
NIST  Standard  Hygrometer.  Second,  an  analysis  was  made  of  all  known  possible 
sources  of  error  and  from  this  analysis,  an  estimate  was  derived  for  the 
accuracy  that  could  be  expected  from  the  generator.  The  results  of  the  two 
methods  were  compared. 

B-3-1  INTERCOMPARISON  TESTS 


A 3x3  Graeco-latin  square  experiment  [34]  was  used  to  test  four  variable 
parameters  of  the  two -pressure  humidity  generator.  The  experiment  was 
designed  to  determine  whether  any  of  the  preselected  levels  of  the  parameters 
could  affect  the  accuracy  of  the  generator.  The  four  parameters  which  were 
tested  in  the  experiment  were  the  presaturator  temperature,  the  final 
saturator  temperature,  the  pressure  of  the  saturator,  and  the  test  gas  (air) 
flow. 

Three  levels  for  each  of  the  four  parameters  were  used  in  the  test. 

The  parameters  were  arranged  in  the  following  form: 
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where 


T = the  final  saturator  temperature  and  T_2o  = -20°C,  = 0®C,  T25  = 

25“C; 

At  = the  difference  in  the  temperature  between  the  presaturator  and  the 
final  saturator,  the  subscript  indicating  the  amount  in  degrees 
Celsius  that  the  former  exceed  the  latter; 

P = the  final  saturator  pressure  and  = 10^  Pa,  P2  =2  x 10^  Pa,  P^  = 

5 X 10^  Pa;  and 

F = the  rate  of  test  air  flow  and  F-j  = .03  m^niin,  F^  = .15  m^/min, 

F^q=  .3  mVniin  (at  25*’C  and  10^  Pa). 

Each  entry  of  the  Graeco -latin  square  represents  a run  and  the  four  parameters 
were  maintained  at  the  designated  levels.  The  mixing  ratio  of  the  moist  air 
produced  by  the  two -pressure  generator  was  calculated  by  using  eq.  (45)  for 
each  of  the  runs  and  the  results  were  compared  with  the  value  of  the  mixing 
ratio  as  measured  by  the  NIST  Standard  Hygrometer.  The  NBS  Standard 
Hygrometer  has  a maximum  uncertainty  of  0.13  percent. 

B-3-1-1  RESULTS 

A percentage  difference,  d,  was  obtained  by  using  the  expression  d = 
(r^-rg)/r^  x 100  where  r^  is  the  computed  mixing  ratio  for  the  generator  and  rg 
is  the  mixing  ratio  measured  by  the  NIST  standard  hygrometer.  The  results  of 
the  tests  for  the  3x3  Graeco-latin  square  experiment  are  given  in  table  22. 
Two  values  in  a box  represent  a repeat  run. 

In  the  first  column  of  table  22,  T.2q  is  the  value  of  the  parameter,  T, 
common  to  all  three  boxes  in  that  column  while  the  other  three  parameters  are 
each  represented  at  the  three  different  levels  Indicated.  Similarly,  T^^  and 
T25  are  the  values  of  the  parameter,  T common  to  all  boxes  in  columns  two  and 
three,  respectively.  If  the  four  tested  parameters  are  assumed  to  be 
independent  of  each  other  and  have  no  interactions,  then  the  average  value  in 
each  column  is  indicative  of  the  correlation  of  the  percentage  difference,  d, 
with  the  value  of  T for  that  column.  Similar  analyses  of  the  rows  yield  the 
correlation  of  d with  F,  analyses  of  one  set  of  diagonals  give  the  correlation 
of  d with  P2,  while  analyses  of  the  second  set  of  diagonals  give  the 
correlation  of  d with  At.  These  results  are  given  in  table  23. 

In  order  to  assess  the  significance  of  d in  table  23,  it  is  necessary  to 
remember  that  the  maximum  uncertainty  in  r^  is  0.13  percent.  Therefore,  if  d 
exceeds  0.13  percent,  the  increase  is  ascribed  to  r^  and,  more  particularly, 
to  the  corresponding  parameter  level. 

B-3-2  ERROR  ANALYSIS 

An  estimate  of  the  maximum  uncertainty  in  the  calculated  mixing  ratio, 
r^,  of  the  generator  for  each  of  the  runs  was  obtained  by  using  the  estimated 
systematic  uncertainty  plus  three  times  the  random  uncertainty  of  the  measured 
pressure  and  temperature  in  the  final  saturator.  The  maximum  uncertainty  of 
the  enhancement  factor,  which  also  includes  the  uncertainty  in  the  saturation 
vapor  pressure  values,  was  obtained  from  table  9 of  Hyland's  paper  [31].  The 
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standard  deviations  were  computed  for  the  measured  temperature  and  pressure  of 
the  saturator  in  each  of  the  runs  and  the  standard  deviation  of  the  mean  was 
used  as  a measure  of  the  random  uncertainty.  The  systematic  uncertainty  in  the 
pressure  measurement  was  attributed  to  the  uncertainty  in  the  calibration  of 
the  pressure  gage,  and  for  the  temperature  measurement,  the  systematic 
uncertainty  was  primarily  due  to  the  maximum  temperature  gradient  detected  in 
the  final  bath,  which  contains  the  saturator  and  the  test  chamber.  The 
estimated  uncertainty  in  the  calculated  mixing  ratio  was  obtained  from  the 
expression: 


P -fe 
s w 


Ae„  AP,  Af 


(51) 


Table  24  lists  the  estimated  systematic  uncertainty  and  the  3o  errors 
for  each  of  the  runs  and  also  the  estimated  maximum  uncertainty  in  the 
generator . 

By  inserting  the  values  of  the  estimated  maximum  uncertainty  given  in 
table  24  in  the  appropriate  squares  of  the  3x3  Graeco -latin  square  and  by 
computing  the  mean  for  each  row,  column  and  diagonal  of  the  square,  the 
results  obtained  are  a measure  of  the  estimated  maximum  uncertainty  of  the 
generator  for  each  level  of  the  tested  parameters.  These  results  are  given  in 
table  25. 

Comparison  of  the  results  (to  the  nearest  tenth  of  a percent)  of  tables 
23  and  25  show  that  the  calculated  maximum  uncertainty  of  the  generator  is 
equal  to  or  greater  than  the  percentage  difference  d for  the  various  levels  of 
the  tested  parameters  and  therefore  it  may  be  concluded  that  the  performance 
of  the  generator  is  not  affected  by  the  four  tested  parameters  over  the  range 
in  which  these  parameters  were  tested. 

B-4  SUMMARY  OF  UNCERTAINTY  QUOTE  FOR  THE  NIST  TWO-PRESSURE  HUMIDITY 

GENERATOR.  MARK  2 

The  results  of  the  intercomparison  tests  of  the  generator  with  the 
standard  hygrometer  which  are  shown  in  the  previous  section  indicate  that  the 
estimated  maximum  uncertainty  of  the  generator  (based  on  the  systematic  and  3o 
uncertainties  of  the  temperature  and  pressure  measurements  and  the  maximum 
error  for  the  enhancement  factor)  are  equal  to  or  greater  than  the  measured 
difference,  d.  Therefore,  similar  calculations  were  made  beyond  the  range 
covered  by  the  intercomparison  tests  to  obtain  the  estimates  of  the  maximum 
uncertainty  for  the  generator  over  the  temperature  range  of  -55  to  80 °C  and 
for  pressures  from  ambient  to  3.3  x 10^  Pa.  Table  26  lists  these  results 
which  are  given  in  units  of  mixing  ratio,  volume  ratio,  dew-point  temperature, 
and  relative  humidity. 

B-5  SURVEILLANCE 
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As  stated  in  Section  B-2,  the  calculations  of  the  various  units  of 
humidity  require  the  measurement  of  the  temperatures  and  pressures  in  the 
final  saturator  and  the  test  chamber. 

To  assure  that  the  accuracy  of  the  temperature  measurements  are 
maintained,  the  (ice  point  resistances)  of  the  platinum  resistance 
thermometers  are  redetermined  at  six  month  intervals  through  use  of  a triple 
point  of  water  cell.  If  the  thermometer  reading  deviates  a few  millidegrees 
from  the  calibrated  value,  the  thermometer  is  returned  to  the  Temperature 
Standards  Group  for  recalibration. 

It  has  been  determined  that  by  checking  the  zero  of  the  pressure  gages 
each  month,  the  required  accuracy  for  the  pressure  measurements  is  maintained. 
In  addition,  the  gages  are  checked  once  a year  against  a calibrated  dead 
weight  piston  gage. 

In  addition,  tests  are  continuously  made  to  check  the  overall 
performance  of  the  generator.  An  example  of  such  a test  is  to  maintain  a 
constant  temperature  and  pressure  in  the  final  saturator  and  to  vary  the 
external  saturator  temperature.  If  a hygrometer  with  high  sensitivity 
indicates  no  change  in  the  output,  the  test  indicates  satisfactory  operation 
of  the  saturation  system.  Similar  tests  are  made  to  determine  the 
independence  of  the  generator  to  air  flow  rates  and  presaturator  water  level. 

Although  there  were  no  indications  regarding  the  malfunctioning  of  the 
generator  after  a period  of  ten  years  since  the  original  test,  an 
intercomparison  test  of  the  generator  was  recently  made  with  the  NIST  Standard 
Hygrometer.  The  test  was  performed  at  a mixing  ratio  of  10  mg/g  and  the 
predicted  mixing  ratio  of  the  generator  agreed  to  within  .01  percent  of  the 
measured  mixing  ratio  of  the  gravimetric  hygrometer.  The  result  of  the 
intercomparison  test  reconfirmed  that  the  overall  operation  of  the  generator 
has  not  changed. 
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Table  1.  Standard  deviation  of  the  air  density  for  water 
vapor  mass  determination. 


Ind. 

Partial 

S tandard 

. Deviation 

Variance 

Var . 

Derivative 

dz 

s(Zn) 

1 az  / ^ 

/ ^ 

)s2(Zn) 

n 

\ V 

V "y 

1 

T 

3.9  X 10"^g/cm^/K 

0.07K 

0.27  X 10*^g/cm^ 

0.0729 

X 10'''^g^/cm‘^ 

B 

1.6  X lO’^g/cm^/mmHg 

0 . 09  mmHg 

0.14  X lO'^g/cm^ 

0.0196 

X 10""'^g^/cm^ 

RH 

0.15  X 10'6g/cm2/%RH 

0.77oRH 

0.10  X lO'^g/cm^ 

0.0100 

X 10'*'^g^/cm*^ 

0 . 24  X lO'^g/cm^/mmHg 

0.09  mmHg 

0.022  X 10‘^g/cm^ 

0.0005 

X lO'^^g^/cm*^ 

o 

II 

Q. 

(M 

w 

,1030  X 

10*^2g2/cm6 

s(p)  = 0.32  X 10*^g/cm^ 
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Table  2.  Standard  deviation  of  external  u-tube  volume. 


Ind.  Partial  Standard  Deviation  Variance 

Var.  Derivative 


3^0 

S(z„) 

Pi 

53 . 9cm^/g 

0.42  X lO'^g/cm- 

P2 

3 . 8cm^/g 

0.42  X lO’^g/cm- 

2 . 6cm^/g 

45  X 10‘^g/cm' 

68 . Icm^/g 

45  X lO'^g/cm- 

Wl 

1 . Ocm^/g 

0.013g 

“2 

1 . Ocm^/g 

O.OlOg 

U) 

22.6  X lO'^cm^ 

negligible 

1.7  X lO'^cm^ 

negligible 

117  X lO'^cra^ 

negligible 

3064  X lO’^cm^ 

9.4  X 10‘^cm^ 

0.013cm2 

169  X lO’^cm^ 

0 . OlOcm^ 

100  X lO'^cm^ 

= 278  X lO'^cm^ 

s(V„) 

= 0.017cm^ 
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Table  3.  Standard  deviation  of  the  buoyancy  correction 
for  the  water  vapor  mass  determination. 


Ind . 
Var. 

Partial 

Derivative 

Standard 

Deviation 

Variance 

3% 

f 

Pi 

11.36cm2 

0.32  X 10‘^g/cm^ 

1,  V 

3.635  X lO'^g 

\^W 

13.213  X 10'''2g2 

Pf 

11 . 36cm^ 

0.32  X 10*^g/cm^ 

3.635  X lO'^g 

13.213  X 10*^2g2 

''t 

lO'^g/cm^ 

0 . 04cm^ 

4.00  X lO'^g 

16.000  X 10b*''2g2 

lO'^g/cm^ 

0 . 04cm^ 

4.00  X 10‘6g 

16.000  X 10'‘'2g2 

«Li 

0.00015 

5.2  X lO'^g 

Negligible 

Negligible 

“if 

0.00015 

5.2  X 10'6g 

Negligible 

Negligible 

w 

^Si 

0.00014 

lO-'^g 

Negligible 

Negligible 

W 

^Sf 

0.00014 

lO-'^g 

Negligible 

Negligible 

s2(C^)  = 58.426  X 10'‘'2g2 
s(C^)  = 7.6  X lO’^g 
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Table  4.  Standard  deviation  of  hydrogen  density  in  sealed  U-tube. 


Ind.  Partial  Standard  Deviation  Variance 

Var.  Derivative 


SPh 

S(z^) 

T 

0.27 

X lO'^g/cm^/K 

0.07K 

1.9  X 10*® 

3.6  X 10**’®g^/cm® 

B 

0.11 

X lO’^g/cm^/minHg 

0 . OSminHg 

0.66  X 10*® 

0.44  X 10*'’®g^/cm® 

s^(p^)  ^ 4.0  X lO'^^g^/cm^ 
s(p^)  = 2 X 10’®g/cm^ 
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Table  5.  Standard  deviation  of  initial  U-tube  volume. 


Ind.  Partial  Standard  Deviation  Variance 

Var.  Derivative 


Wg  0.93  X lO^cm^/g 

0.93  X lO^cm^/g 
Pg  2.9  X lO^cm^/g 

2.9  X lO^cm^/g 


80  X lO'^g 

80  X 10*6g 

0.29  X lO'^g/cm^ 
0.02  X lO'^g/cm^ 


74 

X 

10‘^cm^ 

74 

X 

lO'^cm^ 

8 

X 

lO'^cm^ 

0.6 

X 

10‘^cm^ 

/avA^ 


dz. 


s"(Zn) 


5476  X lO'^cm^ 
5476  X lO'^cm^ 
64  X lO'^cm^ 


Negligible 


s2(V.)  = 11016  X lO'^cm^ 
s(V.)  = O.llcm^ 
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Table  6.  Standard  deviation  of  the  correction  for  mass  of  internal  gas. 


Ind.  Partial  Standard  Deviation  Variance 

Var.  Derivative 


ac, 

s(2n) 

1 

u:)‘" 

1 X 10‘^g/cm^ 

0 . llcm^ 

1.1 

X lO'^g 

1.2  X 10*''°g2 

V* 

1.17  X lO'^g/cm^ 

0 . 066cm^ 

7.7 

X 10‘^g 

59.3  X lO-'lOg^ 

Pf 

30cm^ 

0.29  X lO'^g/cm^ 

9 

X lO'^g 

0.8  X lO-'i^g^ 

Pi 

29.5  cm^ 

0.29  X lO'^g/cm^ 

9 

X lO’^g 

0.8  X 10’'’°g2 

s^(C  ) = 62.1  X 10''°g2 
s(Cg)  = 7.9  X 10-5g 


’'^In  the  case  of  the  second  and  third  tubes , the  error  in  v has 
been  shown  to  be  negligible.  For  these  tubes,  it  can  be  shown 
that  s(C  ) = 1.7  X lO'^g. 
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Table  7.  Standard  deviation  of  the  mass  of  water  collected  in  a u-tube*. 


Independent 

Variable Standard  Deviation Variance 


s(Zn) 

s2(z„) 

Wf 

11  X 10*6g 

121  X 10'^2g2 

w.- 

11  X lO'^g 

121  X 10-''2g2 

7.5  X lO'^g 

56  X 10-^2g2 

Cb 

7.6  X lO'^g 

58  X 10*‘'2g2 

Ca 

Negligible 

Negigible 

c ** 
g 

79  X lO'^g 

6241  X 10‘‘'2g2 

s^(m)  = 

6597  X 10*‘'2g2 

s (m)  = 

81  X 10*6g 

*The  error 

arising  from  the  non- application  of  C-  is 

considered  on  p.  423 

’^■^Note  that  there  is  another  case,  where  s(Cg)  is  1.7  x lO'^g  (see  Table  6). 
In  that  case,  it  can  be  shown  that  s(m)  = 25  x lO'^g. 
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Table  8. 

Dry 
of  a 

runs  - estimate 
single  tube. 

of  the  accuracy  of  the  weight  gain 

Tube 

No. 

Run 

No. 

Average 

U-Tube 

Weight* 

(mg  X 10^) 

Difference 

in  Average  Degrees  of 

Weight  (Difference)^  Freedom 

(mg  X 10^)  (nig^  X 10^) 

11 

19 

1 

35 

+16 

256 

2 

25 

-10 

100 

3 

30 

+ 5 

25 

4 

16 

-14 

196 

5 

17 

+ 1 

1 

6 

17 

0 

0 

5 

12 

19 

7 

23 

+ 4 

16 

8 

29 

+ 6 

36 

9 

34 

+ 5 

25 

10 

28 

- 6 

36 

11 

23 

- 5 

25 

12 

21 

- 2 

4 

5 

14 

13 

13 

14 

+ 1 

1 

14 

8 

- 6 

36 

2 

15 

9 

+ 1 

1 

E (Difference)^  = 758  E = 12 

Variance  = 758  = 632  x 10‘^rag^ 

12 

Standard  Deviation  = 8.0  x lO'^mg 


’'^Only  the  final  significant  figures  are  given  since 
only  these  changed  from  run  to  run. 
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Table  9.  Moisture  pick-up  experiment. 


Run 

No. 

Length 

of 

Run 
(min. ) 

Tube 

No. 

Difference  Between  the 
Change  Mass  Increase  in  the 

in  Drying  Train  and  the 

Mass*  Mass  Decrease  in  the 

Water  Source+++ 

(mg  X 10^) 

(Difference)^ 

(mg^  X lO'*) 

1 

225 

14** 

11+ 

+1.12445 
+ .00107 

+1.12552 

15++ 

-1.12532 

+20 

400 

2 

225 

14** 

11+ 

+1.09609 
+ .01721 

+1.11330 

15++ 

-1.11349 

-19 

361 

3 

300 

14** 

11+ 

+1.15937 
+ .01015 

+1.16952 

13++ 

+1.16950 

+ 2 

4 

4 

296 

J** 

R+ 

+ .74388 
+ .00011 

+ .74399 

A++ 

- .74407 

- 8 

64 

5 

490 

B** 

J+ 

+1.26387 
+ .00008 

+1.26395 

R++ 

-1.26381 

+14 

196 

E 

Variance 
Standard  Deviation 

= 1025mg2  X 10’^ 

= 256mg^  X 10*^ 

= 16mg  X 10'^ 

'A’Based  on  average  of  two  to  five  weighings  before  and  after  the  run. 
**First  tube  in  drying  train. 

-t-Second  tube  in  drying  train. 

-f“i-Water  souce . 

+++Plus  sign  indicates  drying  train  mass  increase  was  greate  than  the  water 
source  mass  decrease. 
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Table  10.  Internal  volumes  of  cylinders  at  250° C. 


Cvlinder 

No.  1 

Cvlinder 

No.  2 

Run  No 

Volume, 

7 ° 

(cm^) 

Deviation 
from  Mean 
(cm^) 

Volume,  V 

7 ° 

(cm^) 

Deviation 
from  Mean 
(cm^) 

2 

29658.8 

-2.2 

29734.2 

-0.8 

3 

29647.6* 

- 

29732.9 

-2.1 

4 

29654.0 

-7.0 

29736.6 

+1.6 

5 

29661.4** 

+0  A 

29738.1 

+3.1 

6 

29665.5 

+U.5 

29733.4 

-1.6 

7 

29665.3 

-^4.3 

- 

- 

Mean 

29661.0 

±3.68 

29735.0 

±1.84 

Standard 

deviation  of  the  mean 

±2.15 

0.99 

Maximum 

-7.0 

+3.1 

*Air  trapped  in  cylinder.  This  value  was  not  used  to  compute  mean. 

'^^Corrected  for  water  trapped  in  vacuum  exhaust  valve. 
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Table  11,  Standard  deviation  of  the  buoyancy  correction  in  the 
cylinder  volume  determination. 


Ind.  Partial  Standard  Deviation  Variance 

Var.  Derivative 


s(Zn) 


s"(Zn) 


Pi 

280cm3 

0.42  X 10*^g/cm^ 

O.OOOg 

0.000000g2 

pf 

26,420cm^ 

0.42  X lO'^g/cm^ 

O.Ollg 

0.000121g2 

PH20 

34 . 5cm^ 

45  X lO'^g/cm^ 

0.002g 

0.000004g2 

Wf 

1.00103 

0.050g 

O.OSOg 

0.002500g2 

1.00964 

O.lOOg 

O.lOlg 

0.010201g2 

s2(c^^)  = 0.012826g2 

s(C^,)  = 0.113g 
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Table  12.  Standard  deviation  of  the  mass  of  water  used 
in  the  cylinder  calibration. 


Independent 

Variable 

Standard 

Deviation 

S(z„) 

Variance 

s2(z„) 

Wf 

0.040g 

o.ooieoog' 

Wf 

0.020g 

0.000400g' 

0.030g 

0.000900g' 

Cr 

0.005g 

0.000025g‘ 

Cb 

0.113g 

0.012770g' 

s2(M)  = 0.015695g2 

s(M)  = 0.125g 
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Table  13.  Standard  deviation  of  the  cylinder  volume  at  the 
calibration  temperature. 


Ind.  Partial  Standard  Deviation  Variance 

Var.  Derivative 


av. 


az. 


s(z^) 


s"(Zn) 


M 1 cm^/g  0.125g 

d 29.6  X lO^g  45  x lO'^g/cm^ 


0 . 125cra^ 
1 . 33cra^ 


0.0157cm^ 

1.77cra^ 


s2(V^)  = 1.79cm^ 
s(V^)  = 1.34cm^ 
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Table  14,  Standard  deviation  of  the  cylinder  volume 
correction  to  a standard  temperature. 


Ind.  Partial  Standard  Deviation  Variance 

Var.  Derivative 


z 


n 


S(z„) 


sHz„) 


V. 


93  X lO'^cm^/cm^ 
59.4  X 10^cmV°C 
1.38cmV°C 


1 . 34cm^ 

0.45  X 10'^cm^/cm^/° C 
0.1°C 


0 . 000cm  4^ 
0.027cm3 
0 . 138cm^ 


0 . OOOOcm^ 

0 . 0007cm^ 

0 . 0190cm^ 


s^(C^)  - 0.0197cm^ 
s(C^)  ^ 0 . 14cm^ 


52 


Table  15.  Standard  deviation  of  p . 

^ 0 


Ind. 

Var. 

Partial 

Derivative 

Standard  Deviation 

Variance 

^Po 

S(z„) 

f 

s2(z„) 

dz_ 

V 

1 az  / 

n 

\ 7 

\ n/ 

%N2 

1.25064  X 10*^g/cm^ 

0.00004 

5.0026  X 10‘®g/cm2 

25.026  X 10-‘'6gVcm^ 

ZOj 

1.42848  X lO’^g/cm^ 

0.00002 

2.8570  X lO-Sg/cm^ 

8.162  X 10*'’^g^/cm^ 

%C02 

1.96465  X lO'^g/cm^ 

0.00004 

7.8586  X 10'®g/cm^ 

61.758  X lO’^^gVcm^ 

%A 

1.78533  X lO'^g/cm^ 

0.00001 

1.7853  X 10'®g/cra2 

3.187  X 10‘^^g^/cm^ 

»a 

2.678  X 10*®cm'2 

6 X lO'^g 

16.0680  X 10'''2g/cm2 

Negligible 

V 

5.769  X lO-Sg/cm^ 

0 . 9cm^ 

5.1921  X 10'®g/cm^ 

26.958  X 10*‘'6g2/cm^ 

s^(Pq)  = 0.0125  X lO'^^g^/cm^ 
s(Pq)  = 0.1  X 10“^g/cni^ 
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Table  16.  Standard  deviation  of  the  density  of  air  in  a cylinder. 


Ind.  Partial  Standard  Deviation  Variance 

Var.  Derivative 


z 

1.09  X lO'^g/cm^ 

54  X 10*^ 

58.9 

X 

lO’^g/cm^ 

0.0035 

X 

lO'^^g/cm^ 

Po 

0.844 

0.1  X 10*^g/cm^ 

0.0844 

X 

lO'^g/cm^ 

0.0071 

X 

10*^^g/cm^ 

T 

3.67  X lO’^g/cm^K 

0.008K 

0.029 

X 

10‘^g/cm^ 

0.0008 

X 

10'*’^g/cm^ 

B 

1.57  X lO'^g/cm^/nuiiHg 

0.06  Torr 

0.09 

X 

lO’^g/cm^ 

0.0081 

X 

10-‘'2g/cm6 

s^(p)  = 0.0195  X lO’^^g^/cm^ 
s(p)  = 0.14  X lO’^g/cm^ 
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Table  17.  Relative  systematic  error  in  mixing  ratio  from 
incomplete  water  absorption*. 


No.  of  Cylinders 

Nominal 

Mixing  Ratio  Relative  Error 

(mg/g)  (parts/10^) 

0.7 

27.1  -0.017 

1 

18.9  -0.025 

10 

1.89  -0.25 

40 

0.475  -0.99 

100 

0.189  -2.5 

*Mass  water  vapor  = 
air  per  cylinder  = 

600mg,  regardless  of  duration;  mass  dry 

31. 6g. 
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Table  18.  Water  vapor  and  air  masses  collected  from 
leakage  and  the  test  gas. 


No . of 
Cylinders 

Room 

Temperature 
and  RH 

A B 

Water  Gained  From 

Leak  Test  Gas 

(mg)  (mg) 

C 

Air 

Leak 

(g) 

D 

Gained  From 

Test  Gas 

(g) 

0.7 

19°C 

6.9  X 10-5 

600 

4.91 

X 

10-5 

22.1 

1 

10%  RH 

9.8  X 10'5 

600 

7.02 

X 

10-5 

31.6 

10 

98  X 10*5 

600 

70.2 

X 

10'5 

316 

40 

392  X 10'5 

600 

280 

X 

10-5 

1264 

100 

980  X 10'5 

600 

702 

X 

10-5 

3160 

0.7 

28°C 

69  X 10-5 

600 

4.91 

X 

10-5 

22.1 

1 

60%  RH 

98  X 10’5 

600 

7.02 

X 

10'5 

31.6 

10 

980  X 10'^ 

600 

70.2 

X 

10'5 

316 

40 

0.039 

600 

280 

X 

10’5 

1264 

100 

0.098 

600 

702 

X 

10'5 

3160 
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Table  19.  Relative  systematic  error  in  mixing  ratio  from  leakage, 
for  runs  of  various  lengths . 


No . of 
Cylinders 

Room 

Temperature 
and  RH 

E 

"True"  Mixing 
Ratio  (mg/g)* 

F 

Observed  Mixing 
Ratio  (mg/g)+ 

Error  (parts/10' 
[(F-E)/E]  X 10^ 

0.7 

19°C 

27.14932 

27.14926 

-0.02 

1 

10%RH 

18.98734 

18.98730 

-0.02 

10 

1.898734 

1.898733 

-0.01 

40 

0.4746835 

0.4746856 

+0.04 

100 

0.1898734 

0.1898761 

+0.14 

0.7 

28°C 

27.14932 

27.14929 

-0.01 

1 

60%RH 

18.98734 

18.98733 

+0.01 

10 

1.898734 

1.898761 

+0 . 14 

40 

0.4746835 

0.4747133 

+0.63 

100 

0.1898734 

0.1899040 

+1.61 

^Column  B/D,  Table  18 


+Columns  (A+B)/(C+D),  Table  18 
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Table  20.  Summary  of  maximum  errors  (parts  per  10^). 


No . of 
Cylinders 

Nominal  Mixing 
Ratio  (mg/g) 

3 X SD 

Errors  in 
M,p,V 

Incomplete  Water 
Absorption 

Leakage 

0.7 

27.1 

±8.1 

±2.09 

-0.017 

-0.02 

1 

18.9 

±8.1 

±2.09 

-0.025 

-0.02 

10 

1.89 

±8.1 

±2.09 

-0.25 

-1-0.14,  -0.01 

40 

0.475 

±8.1 

±2.09 

-0.99 

-f-0.63 

100 

0.189 

±8.1 

±2.09 

-2.5 

-Hi.  61 

Table  21.  Error  band  in  mixing  ratios  . 


No.  of  Cylinders 


Mixing  Ratio  Error  Band 

(mg/g)  (parts/10^) 


0.7 

27.1 

CNI 

o 
!— 1 

to 

+10.2 

1 

18.9 

-10.2 

to 

+10.2 

10 

1.89 

-10.4 

to 

+10.3 

40 

0.475 

-11.2 

to 

+10.8 

100 

0.189 

-12.7 

to 

+11.8 
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Table  22.  The  percentage  difference  between  r^  and  r^ 


Table  23,  The  correlation  of  parameter  levels  with  the  percentage 

difference  between  r and  r 

w g 


d 

T.20  0.23% 

d 

At^5  0.10% 

d 

P5  0.10% 

d 

F,(,  0.17% 

T^  0.11% 

At^  0.18% 

P2  0.18% 

F5  0.17% 

T25  0.12% 

At  0.17% 

0 

P^  0.18% 

0.14% 
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Table  24.  Experimental  Errors 


< 


fA 

<1 


Q) 

<1 
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Table  25. 


The  estimated  maximum  uncertainty  of  the  generator's  mixing 
ratio  for  the  designated  levels  of  the  tested  parameter 


'^-20 

0.33% 

At, 5 

0.19% 

^5 

0.22% 

Tio 

0.20% 

To 

0.14% 

Ats 

0.19% 

^2 

0.15% 

^5 

0.21% 

'^25 

0.10% 

rt 

o 

0.18% 

0.22% 

0.16% 

Table  26.  NBS  Two -Pressure  Humidity  Generator,  Mark  2,  Range  and  Accuracy 


Humidity  parameter 


Range 


Accuracy^ 


Mixing  ratio,  r^  (g  water  vapor/kg 

dry  air) 


Volume  ratio,  V (ppm) 


Dew-point  temperature,  T^  °C 


Relative  humidity,  RH  (%)  at  test 

chamber  temperature  T^  (°C)  of: 

-55  ^ T^  < -40 

-40  ^ T^  < -20 

-20  s T^  < 0 

0 ^ T^  < +80 

c 


0.0005 

0.0015 

3.0% 

of 

value 

0.0015 

0.005 

1.5% 

of 

value 

0.005 

0.1 

1.0% 

of 

value 

0.1 

0.3 

0.5% 

of 

value 

0.3 

^ < 

515 

0.3% 

of 

value 

V 

> 

vt 

3 

3.0% 

of 

value 

3 s V < 

10 

1.5% 

of 

value 

V 

> 

V) 

o 

170 

1.0% 

of 

value 

170  ^ V < 

500 

0.5% 

of 

value 

500  s V < 

820,000 

0.3% 

of 

value 

-80  ^ T^ 

< -70 

0.2 

-70  s T^ 

< -35 

0.1 

-35  ^ T, 

< +80 

0.04 

3-98 

1, 

.5 

3-98 

0, 

.8 

3-98 

0, 

.4 

3-98 

0, 

,2 

^The  estimated  bounds  to  systematic  error  plus  three  times  the  standard 
deviation. 
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